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Abstract: This work considers the way that quantum loop effects modify the prop- 
agation of hght in curved space. The calculation of the refractive index for scalar 
QED is reviewed and then extended for the first time to QED with spinor particles 
in the loop. It is shown how, in both cases, the low frequency phase velocity can 
be greater than c, as found originally by Drummond and Hathrell, but causality is 
respected in the sense that retarded Green functions vanish outside the lightcone. 
A "phenomenology" of the refractive index is then presented for black holes, FRW 
universes and gravitational waves. In some cases, some of the polarization states 
propagate with a refractive index having a negative imaginary part indicating a po- 
tential breakdown of the optical theorem in curved space and possible instabilities. 



1. Introduction 

The quantum theory of photon propagation in curved spacetime raises many chal- 
lenging questions, particularly in relation to the realization of causality and unitarity 
in quantum field theories with a fixed gravitational background. In particular, in pre- 
vious work [1-3] we have shown that vacuum polarization effects in curved spacetime 
lead to the phenomenon that curved spacetime acts like an optical medium and can 
be described by means of a refractive index. However, this index has a novel ana- 
lytic structure, compared with conventional optical media, and standard formulae, 
notably the Kramers- Kronig dispersion relation are violated. This analytic structure 
has an origin in the geometry of null geodesic congruences and their conjugate points 
which give rise to curved-spacetime-specific singularities in the Van Vleck-Morette 
determinant and consequently the Green functions. 

This line of research began with the apparent paradox of how to reconcile the 
requirements of causality with the observation by Drummond and Hathrell [4] that 
at low frequencies the phase velocity of light in a gravitational background can be 
superluminal {i.e. the refractive index n{uj) < 1 for uj small). Since causality plays 
such an important role in the theory, it is worthwhile stating exactly what is meant in 
the context of a QFT. Causality in this context is often called "micro-causality", and 
is the property that the retarded/advanced Green functions only have non- vanishing 
support inside (or on) the backward/forward lightcone. The requirement of causality 
can be linked directly to the refractive index itself [1] leading to two conditions: (i) 
n{u!) 1 for large |a;|; and (ii) n(a;) must be an analytic function in the upper-half 
of the complex u plane. The first condition is equivalent to the requirement that 
the wavefront velocity, the high frequency limit of the phase velocity, is equal to 
c [5,6]. Notice that neither of these conditions could be tested in the original work 
of Drummond and Hathrell which was restricted to small uj. However, prior to the 
unlocking of the full frequency-dependence of n{uj) in [1-3], the following argument 
led to a paradox. In conventional dispersive, or dielectric, media the refractive index 
obeys the familiar Kramers-Kronig dispersion relation 

n(oo) = n(0) — [ — Im n{uj) . (1.1) 
Jo ^ 

Since Imn(a;) is constrained to be positive, being related by the optical theorem (itself 
a consequence of unitarity) to a forward scattering cross-section, eq.(|l.l|) implies 
n{oo) < n(0). This would mean a superluminal low frequency refractive index ?t,(0) < 
1 would be incompatible with a causal limit n{oo) = 1. The resolution, found in 
refs. [1-3], is that spacetime behaves like an optical medium for which the Kramers- 
Kronig relation no longer holds in the form ( p.. ID and, furthermore, Im?T,(a;) is not 
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constrained to be positive. One can trace this non-compliance to the novel analytic 
structure — compared with conventional dielectric media — of the refractive index of 
curved spacetime. 



In curved spacetime, the refractive index will in general be a local quantity 
n{x,uj), since spacetime is not generally homogeneous, (as well as a matrix quantity 
in order to describe the two polarization states of the photon). What we will find is 
that n{x, uj) is consistent with micro-causality , as a consequence of analyticity in the 
upper-half uj plane, and consequently one can write a version of the Kramers- Kronig 
relation in the form 



which does hold. The usual relations that hold in models of dispersive media, namely 
the fact that n{uj) is an even function of u and the property of hermitian analyticity, 
n{uj — ie) = n{uj + ze)*, which together imply that n{(jj + ie) = n{—uj + ie)*, then 



see, allows the causality constraint n(x, oo) = 1 to hold even in the presence of a 
superluminal n(a;, 0) < 1. 

Although this resolution of causality with a low-frequency superluminal phase 
velocity is achieved through the modification of the Kramers-Kronig relation to ac- 
commodate the novel analytic structure of n{x, uj) in curved spacetime, in the course 
of studying particular examples in refs. [1-3] we found cases where Imn(x, tu) is in 
fact negative. This is unexpected since one would have expected this imaginary part 
to signal the decay of the photon into real e^e~ pairs and hence be strictly posi- 
tive: it appears that the optical theorem cannot hold — at least in its conventional 
form — in curved spacetime [7]. The correct interpretation of a negative imaginary 
part, the fate of the optical theorem and implication that the photon modes have an 
increasing amplitude, will be left to future work; one of the purposes of the present 
paper is to collect the "phenomenological" evidence on this issue by studying several 
relevant spacetimes. 

In order to establish these results, we calculated the complete frequency depen- 
dence of the refractive index for scalar QED, obtaining a concise formula involving 
the Van Vleck-Morette (VVM) matrix [1].^ This makes clear the geometric origin of 
the final result, and therefore its generality: loop corrections to Green functions in 
curved spacetime will generically have a richer analytic structure than in flat space- 
time. The quantum corrections to photon propagation are governed by the one-loop 

Usually one speaks of the Van Vleck-Morette determinant; here, we are interested in the actual 
matrix which we also call the "VVM matrix" . 




(1.2) 



lead to ( |1 . 1| ) . Neither of these relations hold in curved spacetime which, as we shall 
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vacuum polarization. As shown in rcfs. [1-3], both in a worldhne formahsm and 
using the more conventional heat kernel (or proper time) formalism, when the mass 
of the electron is much greater than the curvature scale and so the "geodesic ap- 
proximation" applies, this is determined by a path integral over fluctuations around 
the classical null geodesic traced by the photon; in turn, these are governed by the 
geometry of geodesic deviation, explaining the origin of the VVM matrix in the final 
result. The great simplification is that to leading order in a weak curvature expan- 
sion, the geometry of geodesic deviation is entirely encoded in the Penrose limit [8] of 
the original spacetime around the photon's classical null trajectory. It follows that 
the whole problem of determining the refractive index can be resolved simply by 
studying the appropriate Penrose plane wave limit, with all the relevant information 
on the background spacetime being encoded in the plane wave profile function hij{u) 
to be defined in due course. 

In this paper, we develop these fundamental ideas in a number of directions, in 
effect building up a phenomenology of the interesting and counter-intuitive effects 
that can arise in what we have previously termed "quantum gravitational optics" [9]. 
We begin by reviewing and simplifying the conceptual basis of our calculation of the 
refractive index for scalar QED in curved spacetime, highlighting the role played 
by the Penrose limit in determining the vacuum polarization and the relation of 
the geometry of null geodesic congruences to the analytic structure of the refractive 
index. Next we calculate the refractive index for spinor QED, introducing the nec- 
essary geometric tools, in particular the spinor parallel transporter in a plane wave 
spacetime, and evaluating the vacuum polarization for a massive spinor loop. The 
result is qualitatively similar to the simpler case of scalar QED and displays the same 
generic analyticity properties, but also shows some physically interesting differences 
in particular models. 

We then consider photon propagation in a variety of examples of gravitational 
backgrounds in order to build up intuition about the role of symmetries, singularities 
and time-dependence of the curved spacetime in the realization of causality and 
unitarity in QED. In each case, we need to find the Penrose limit corresponding to the 
classical photon trajectory of interest, then solve the geodesic deviation equations to 
construct the VVM matrix from which the vacuum polarization and refractive index 
follow. In most cases, we have to resort to a numerical evaluation of the refractive 
index itself. 

The first examples are homogeneous plane waves. These are interesting both as 
simple toy models and as the Penrose limits of certain geodesies in classical black hole 
and cosmological spacetimes near the singularity. They have a large spacetime isom- 
etry group, an extension of the Heisenberg algebra, and can be usefully realized as 
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coset spaces. We consider two cases, the symmetric plane waves originally studied in 
refs. [1-3] and the singular homogeneous plane waves which arise as near-singularity 
Penrose limits. The next set of backgrounds that are considered are the black hole 
spacetimes, especially the Schwarzschild and Kerr metrics. These are examples of 
Petrov type D spacetimes and we can exploit the resulting simplifications to give 
a very general description of the Penrose limits for various choices of geodesies. In 
particular, the Penrose limit corresponding to the principal null geodesies are flat, 
implying that the quantum corrections to the refractive index vanish identically, a 
result already observed in the low- frequency limit. We also see how in general the 
Penrose plane wave profile depends on the Walker-Penrose integral of motion [10,11] 
characterizing the classical trajectory and show how the near-singularity limits re- 
duce to singular homogeneous plane waves in the Penrose limit. 

We then go on to consider FRW universes. In these spaces the definition of 
the refractive index needs a slight modification because the singularity in the past 
does not allow one to consider waves coming in from past null infinity. Finally, we 
discuss the behaviour of the refractive index in both a weak gravitational wave and 
a gravitational Shockwave. 



2. Vacuum polarization and Photon Propagation in Curved 
Spacetime 

We begin by reviewing briefly the eikonal formalism and the derivation of the refrac- 
tive index in the low-frequency effective action from the QED effective action [1] . 

In classical electrodynamics, the propagation of photons in curved spacetime 
is governed by the Maxwell equation V ,j,F^^ = 0, with F^,^ = V^Aj, — Vi^Aj^. In 
the eikonal approximation, the electromagnetic field is written in terms of a slowly- 
varying (with respect to the curvature scale) amplitude A{x) and a rapidly-varying 
phase Q{x) as follows: 

A,{x) = ^(x)£V(x)e^®(^) , (2.1) 

where i^{x) is the polarization tensor. The wave-vector is identified as = 9^0 
and we fix the gauge so that the two independent polarizations i = 1,2 satisfy 
the transverse condition k ■ = 0, have vanishing component along k^ and are 
spacelike normalized such that • = 5^. The eikonal expansion is in powers 
of the frequency uj and is valid in the regime uj ^ y/R, with it! a typical curvature 
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scale. ^ The wave vector itself is 0{uj) whereas A and are 0{1), so at leading 
order, the Maxwell equation gives 



= dQ.dQ = . (2.2) 

Since this also implies k ■ Vk^^ = 0, it follows that the integral curves of the vector 
field k^ are null geodesies, which are identified as the classical trajectories of the 
photon. At next-to-leading order, the eikonal equations describe the variation of the 
amplitude and polarization along these null geodesies: 

fc-V£^ = 0, 

2.3 

k-V logA = -^V -k . 

The second of these relates the change in amplitude to the expansion 6 = W Jt^ (with 
k^ = LJ~^k'^), one of the optical scalars in the Raychaudhuri equations. 

In order to study photon propagation in a general spacetime, it is convenient to 
use a set of coordinates that are specially adapted to the vector field k^. These are 
the Penrose coordinates, which are also known as ''adapted coordinates" , {u,V,Y°-), 
a = 1,2. Here, u is the affine parameter along a null geodesic, V is the associated null 
coordinate and Y°- are two orthogonal spacelike coordinates. This choice corresponds 
to the embedding of the preferred null geodesic 7 (with = = 0), representing 
the classical photon trajectory, in a twist-free null congruence labelled by constant 
V, Y"-. The metric g^^ can always be written in terms of these adapted coordinates 
in the form [12]:^ 

ds'^ = -2du dV + C(m, V, Y^)dV'^ + 2Ca{u, V, Y^)dY'' dV + Cab{u, V, Y^)dY'' dY^ . 

(2.4) 

The eikonal phase is then taken to be 6 = —ujV, the Fourier mode appropriate for a 
metric with an isometry characterized by a Killing vector dy-, so that A;^ = (1, 0, 0, 0). 
In terms of these coordinates the transport equations ( pl3|) are 

V.^'^ = , 

2.5 

du\ogA = -\e, 

where the expansion scalar is 

e=\C'''duCab = d^\og^ , (2.6) 
where C"* denotes the inverse of Cab and g = —detg^^. 



^We can think of R as the magnitude of a typical element of the Riemann tensor, with the 
understanding that derivatives of it would count a factor of ^/R. 

^In this paper, we work with a mostly plus signature for the metric. 
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Going beyond the classical theory and including the effects of vacuum polariza- 
tion, the Maxwell equation is replaced by 



V,F% = 4 y ,/^d''x'Uj;l°'^^x,x')A''{x') 



(2.7) 



where n)j"^°°P(x, x') is the vacuum polarization tensor appearing in the one-loop ef- 



fective action 



ri-^°°P = - / ^/^d^x^/^d^x' A''{x)nl:l°°'P{x,x')A''{x') . (2.8 



In general, we find that the wave vector is no longer null and the physical light cone 
defined by k"^ no longer coincides with the geometric null cones. In physical terms, 
the phase velocity of light is no longer necessarily c and may be either sub- or super- 
luminal [4]. It can also become polarization dependent, i.e. display gravitationally 
induced bi-refringence. 

To accommodate this, we modify the eikonal phase as follows: 

e = -uJ{VS^j - ^ij{x; iu)) . (2.9) 

where we now think of the phase as a 2 x 2 matrix with respect to the polarizations 
With this form for the electromagnetic field A^{x), and using the metric ( p.4| ), 
we find 

V,Fj; = 2^2 ^^..<^;^) ^g.^^-..y ^ ^2.10) 

to leading order in the eikonal expansion. Provided ^ij{x) is perturbatively small, 
i.e. 0{q), this corresponds to a refractive index matrix: 

n(x;c.) = l + ^^. (2.11) 

The phase 'j9(x; uo) is then determined from the leading, 0{uj^) piece of the r.h.s. of 
( P^.IOD evaluated with the eikonal ansatz (|2.lD,(pi3|) for A^{x). Putting all this to- 
gether, we find the following compact expression for the refractive index: 

n,,{x-u) = 5,, dx'^M^) ^ il^^ix) ni-i°°P(x,x') e[;.)(x') e-(^-^') . 

(2.12) 

The polarizations which propagate with well-defined phase velocities are the linear 
combinations of the basis giving the eigenstates of n. 

To see how this works in the simplest case, we consider first the low-frequency 
limit of the refractive index, which exhibits superluminal phase velocities. This can 
be found by considering the modifications to the Maxwell equation following from 
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the leading terms in a derivative expansion of the one-loop effective action and was 
the approach taken in the original work of Drummond and Hathrell [4] . The relevant 
terms in the low-frequency effective action to one loop are 

r = [ d'x^\-lzF,,F^'' + dV.F^'^V.F^ 

J \- 4: ^ ^2.13) 

+ aRF^.F'^" + bR^^F'^^F-'x + cR^uXpF"" F"^" ' 

where the coefficients for both scalar and spinor QED can be read off from [13] and 
are given in [1]. We find 

a = 0, b = c= d = (2.14) 

for scalar QED, while 

a , 13ci a , a , 

^=-7^7^ (2-15) 



1447rm^ ' SQOirm'^ ' ^QOirm'^ ' SOirm'^ 

for spinor QED, reproducing the original Drummond-Hathrell effective action [4]. 

2 

Z = 1 + log ^ is the one-loop wave-function renormalization factor. 

To find the refractive index, we first evaluate the vacuum polarization tensor from 
( p.l3|) . Since this effective action is local, Il^~^°°^{x,x') is proportional to 6{x,x') and 
( p.l2| ) simplifies immediately since, for example, the amplitude factors cancel and the 
polarizations are evaluated at the same point. This immediately gives the general 
result for the low-frequency limit of the refractive index [4, 9] 



n^f °™(0) = - 2bR,Jij - ScRuiuj . (2.16) 
where Ruiuj = Ru\up^\i)^^ljy So for scalar QED we find 



n^--(O) = + ^^^j^ [RuJij + 2i?„»j) , (2.17) 



while for spinor QED, 



n^P--(O) = 5,, - [im^J^j - AR^,^,) , (2.18) 



Notice that the opposite sign of the b coefficient means that scalars and spinors 
respond oppositely to the Ricci curvature. Since the null energy condition requires 
Ruu > 0, this means the low-frequency phase velocity is superluminal for spinors in 
a conformally flat background, but subluminal for scalars. 

The generalization of the QED effective action to all orders in derivatives was 
found in ref. [13,14] and can be used to extend the expression ( |2.16|) for the refractive 
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index njj(O) to a full perturbative expansion of nij{u)) in powers of uj. This effective 
action consists of ^^RFF" operators acted on by functions of the Laplacian, viz. 

+ + b"oR,,F^^F\ + ^R^^^.F'^'^F^^^ 

+ i-(^i?V^F''^V.F^ + hR^uVxF^'^VpF^'' 

+ hR^^V^F^'^VFx^ + hR.^V^V^F^pF^'' + ^.R^^^^V^F^'V^F^'^ 

(2.19) 

In this formula, the o^, 6„, are known "form factor" functions of three operators, 
i.e. 

S:,^aA%^). (2.20) 
V / 

where the first entry V^^^^/m^ acts on the first following term (the curvature), etc. 
The form factors can be extracted from the very general background field effective 
action originally computed by Barvinsky et al. [15] and are described in detail in 
ref. [13,14]. 



Extracting the refractive index from ( 2.19 ) to leading order in the eikonal ap- 



proximation involves a number of subtleties, which are discussed in detail in ref. [14]. 
The result is an expression of the form 

( \ X ^ X r^('^i^k-V\Ruu , (2iujk-V\Ruiuj 

where the constant coefficients are replaced by functions of the operator k ■ V, which 
describes the variation of the curvature tensors along the original null geodesic 7. 
Since B{x) and C{x) are real functions it follows immediately that even perturba- 
tively at 0{uj), a non-constant curvature along 7 will give rise to an imaginary part 
of the refractive index, which can be positive or negative depending on the variation 
of the curvature along the geodesic. 



The expression above ( |2.21| ), captures all the terms which are linear in the cur- 



vature. In the present work, we will calculate an expression for the refractive index 
whose scope is much wider because its sums up all powers of the curvature (as well 
as its derivatives). The only assumptions are that the curvature is weak, in the sense 
that R <^ m^, and the eikonal approximation is valid lu ^ ^/R.^ Schematically we 



*Both these conditions can be dropped if the spacetime is a plane wave. 
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find 

= ^ij + -^Fij ^ (2-22) 



where it! is a generic curvature scale 



5 



To find the fuU frequency dependence of the refractive index, therefore, we need 
to find the exact form of the vacuum polarization n^"j°°P acting on the eikonal ansatz 
for the electromagnetic field. Even the complete expansion of the effective action to 
all-orders in derivatives is not sufficient to capture the non-perturbative behaviour 
of nij{uj) as a function of frequency [2,3,14]. This was achieved in refs. [2,3], in a 
calculation based on the worldline approach to QFT, and subsequently in [1] using 
the more conventional heat kernel, or proper time, representation of the propagators. 

The key insight in both these methods is that when R <^ we can use the 
"geodesic approximation" for the the propagators of the electron and positron in the 
loop. When the loop is coupled to an external photon the geodesic approximation 
leads to a simple picture: the electron and positron follow the original geodesic of the 
photon before annihilating back into the photon. Therefore the vacuum polarization 
is determined by the geometry of geodesic fluctuations about the photon's classical 
trajectory and hence to leading order in "weak" curvature, i.e. 0{R/m'^), the refrac- 
tive index is governed entirely by the Penrose limit of the original curved spacetime. 
This is because the Penrose limit is a truncation of the original spacetime metric 
which captures the tidal forces corresponding to geodesic deviation. Since the Pen- 
rose limit describes a plane wave metric, we find the remarkable simplification that 
at weak curvature, the refractive index for any given spacetime may be calculated 
simply by considering propagation in the associated plane wave background. 

We explain in detail how this comes about in Section 4, where we re-cap the 
essential points of the derivation of the refractive index for scalar QED [1]. First, 
we discuss the most important aspects of the Penrose limit, plane wave spacetimes, 
geodesic deviation and the Van Vleck-Morette determinant from the perspective of 
photon propagation in curved spacetime. 



3. Penrose limit and Geometry of Plane Wave Spacetimes 

In this section, we collect the essential results about the geometry of the Penrose 
^This counting includes powers and derivatives of the Riemann tensor, e.g. VvRvapX counts as 
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limit [8] and plane waves which will be used in the QFT calculations later in the 
paper. For further details and discussion, see the reviews in refs. [1] and [12]. 



3.1 Geodesic deviation and the optical tensors 



Before introducing the Penrose limit, we describe the basic geometry of geodesic 
deviation in the context of the general metric ( |2.4| ) in coordinates adapted to the 
null congruence around 7. We therefore consider the connecting vector which 
connects corresponding points on neighbouring geodesies in the congruence and study 
its evolution along 7. This is determined by the requirement that the Lie derivative 
of along 7 vanishes, i.e. 

Cj^z^" = A; ■ V^^ - (V^F)z^ = 0. (3.1) 

This implies 

VuZ^" = n^z" , (3.2) 

where VL^y = V ^ky, which plays the role of a connection for the Lie derivative, is 
symmetric since the vector field k'^ is a gradient flow. 

To describe geodesic flow and the Raychaudhuri equations, it is sufficient^ to 
consider connecting vectors with transverse components z"" only. In the Penrose 
coordinates (p.4|) , we have'' 

Vuz"" = n\z' , (3.3) 



with 



^ab — -^duCab , (3.4) 



It follows from ( p.3|) that the transverse components of the connecting vector satisfy 
the geodesic deviation equation: 

V^V^^'^ = - R\buz'' = - h^z" , (3.5) 



^Here, we follow the approach presented by Wald [16]. The 4-dim space V of tangent vectors 
is first restricted to a 3-dim subspace V by the condition k ■ z — 0. The final 2-dim vector space V 
is identified as the vector space of equivalence classes of vectors in V with vectors differing only by 
the addition of a multiple of deemed equivalent. In terms of the Penrose (adapted) coordinate 



system (2.4), it is clear that V is realized by restricting to vectors z'^ whose only non-vanishing 
components are z". Similarly for the fundamental tensor field ^ab- Note that Q^^, automatically 
satisfies k^^Vl^^ = by virtue of the null geodesic equation k ■ S/k" — 0. 

''AH the covariant equations in this sub-section are valid in an arbitrary coordinate system. In 
Penrose coordinates, (^^) is simply — 0. In particular, this confirms that itself is a suitable 
connecting vector, with the null geodesies given simply by y" — constant. 
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where^ 

h\ = - v„^]^ - , (3.6) 

that is, 

Raubu = — 2^"^"'' ~'~ '^^uCac C'^'^ duCdb ■ (3.7) 

A solution of ( p. 5] ) is known as a "Jacobi field" on the null geodesic 7. 

The important point here is that the transverse geodesic fluctuations are con- 
trolled entirely by the metric components Cab in (|2.4|) . As we see below, this underlies 
the key role to be played by the Penrose limit. 

The nature of the geodesic flow is most elegantly summarized in the Raychaud- 
huri equations for the optical tensors. These are defined from flab as: 

^ab = -dCab + d'ab + UJab , (3.8) 

where 9, dab, '^ab are respectively the expansion, shear and twist of the null congru- 
ence. Here, the twist is vanishing by construction, since is by definition a gradient 
field, and for this reason (p^) embeds the geodesic 7 in a twist-free null congruence. 

The Raychaudhuri equations are equivalent to ( p.6| ), which reads 

Vul^afe = — flacC'^'^ftdb — Raubu ■ (3.9) 

In terms of the optical scalars: 

dj 

where the Weyl tensor Caubu is the trace-free part of the Riemann tensor ( |3.7| ). 

Notice that the equation for the expansion immediately shows that duO < 
everywhere along the geodesic provided the null energy condition Ruu > holds, 
implying a generic focusing of the null congruence. Moreover, if we define the shear 
scalar as o" = \/2aabd'"'^ and introduce Newman-Penrose scalars $00 = \Ruu and (in 
a null basis m^,fh^ for the transverse directions aligned with the eigenvalues of aab) 
^0 = Cumum, we can write the Raychaudhuri equations as 

duie±a) = -i(^±a)2-2($ooT^o) . (3.11) 

Physically, this shows that at a given point along the geodesic 7, the congruence must 
be focusing in at least one of the transverse directions. While locally, focus/focus 

^In Penrose coordinates, /i°b — ■~du^°'b — ^''c^'^b , or equivalently, hat = — 9„f2ah + flac^'^b ■ 



(3.10) 

9^ab ^aubu 
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and focus/ defocus are permitted, the null energy condition prohibits a congruence 
with defocus/defocus. 

This illustrates a crucial theorem [16], that for a spacetime satisfying the "null 
generic condition" ($00 7^ or ^0 7^ at some point) and the null energy condition, 
every complete null geodesic possesses a pair of conjugate points. Two points p and q 
on a null geodesic are said to be conjugate if there exists a Jacobi field 2;" G V which is 
not identically zero but which vanishes at both p and q. Loosely speaking, this implies 
there is an "infinitcsimally deformed" null geodesic intersecting 7 at both p and q. 
As we see in the next section, this is precisely the property of geodesic fluctuations 
which is reflected in the quantum field theory. Notice though that this definition of 
conjugate points is linked to infinitesimal deviations and does not necessarily lift to 
geodesies of the full metric; in particular, the existence of conjugate points p and q 
does not imply that there is an actual geodesic, other than 7, joining p and q. 



Another important set of coordinates, specifically designed to study fluctuations 
around a given null geodesic, are null Fermi coordinates. Fermi normal coordinates 
are analogues of the familiar Riemann normal coordinates, in which the connec- 
tion coefficients vanish locally, extended from a given point to the whole of a given 
geodesic. The formalism for constructing Fermi coordinates appropriate for a null 
geodesic was developed in ref. [17] and we refer to this paper for further details. 

To construct a null Fermi coordinate system around 7, we begin by introducing 
a pseudo-orthonormal frame = E'^n dx^ {A — u,v,i — 1, 2) such that 

ds^\^ ^ -2E''E'' + SijE'E^ , (3.12) 

which is parallel-transported along 7 and where Eu'^ is the tangent vector. The Fermi 
coordinates = {u,v,y'^) are then essentially the coordinates along the directions 
defined by this frame, with u playing the role of affine parameter. Precisely, we define 

x^ = {u, x'^ = s£;%(wo)x''(0)) , (a = -y, i = 1, 2) , (3.13) 

where x^{s) are geodesies emanating from a point on 7; equivalently, 

dx^ 



dx>^ 



= . (3.14) 



Just as with Riemann normal coordinates, the key property of Fermi coordinates 
is that to linear order in an expansion in powers of the "transverse" coordinates 
x'* around 7, the connection coefficients vanish, r^(^| = (!?((x")^). This allows an 
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expansion of the metric in terms of the Riemann tensor, 

ds^ = — 2dudv + Sijdy^dy^ 

- [RaupuX^x^du^ + ^R^^f^sx'^x^dudx^ + \R^^p5x''x^dxUx^] + 0{{x''f) . 

(3.15) 



Now return to the geodesic deviation equations ( p.3|) and (|3.5|) for the trans- 
verse Jacobi fields, in Fermi coordinates. These take essentially the same form 
except that because of the defining property of Fermi coordinates we can replace the 
covariant derivatives along the geodesic by simple derivatives, i.e. 

-^z' = n\z^ , (3.16) 
du 

with Vtij = Vikj and h^j = R^uju = —du^^j — ^l^k^'^j on 7. 

For deviations around the preferred geodesic 7, the solutions can be written 
in terms of the initial data at some point u' on 7 by introducing matrix functions 
A{u,u') and B{u,u') as follows: 

z\u) = B'j{u,u')z^{u') + A'j{u,u')z^{u'), (3.18) 

where A{u,u') and B{u,u') themselves satisfy the Jacobi equation 

A + hA = 0, B + hB = 0, (3.19) 

with boundary conditions A{u',u') = 0, duA{u,u')\u=u' = 1, B{u',u') = 1, and 
duB{u,u')\u=u' = 0. Clearly, the functions A and B carry the same information 
about the geodesic flow as the optical tensors. We can make this explicit by com- 
bining (IJ) and i ^JEi ) to give the relation: 



B{u,u') + A{u,u')n{u') = exp 
defining the optical tensors as before from Qij. 



du" (^\95ij + dij) 



(3.20) 



For our purposes, an interesting case is to define a congruence, and therefore 
the optical tensors, by choosing "geodesic spray" boundary conditions such that 
z*(m') = 0. This isolates the function A{u,u') and we find 

du\og A{u,u') = n{u) . (3.21) 
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Taking the trace gives the identity 

dulog det A{u,u) = 9{u) . (3.22) 



The importance of this function A{u,u') lies in its relation to the Van-Vleck 
Morette matrix and geodesic interval in a plane wave spacetime. Indeed, solving 
( |3.19|) for A with appropriate boundary conditions proves to be the quickest way to 



determine the VVM matrix and the propagators in the Penrose plane wave limit. 



3.2 Penrose limit and plane wave geometry 



The Penrose limit associates with any metric (^^i, and null geodesic 7 a plane wave 
spacetime which encodes the geometry of geodesic deviation around 7. Formally, it 
is found by making a global Weyl transformation followed by an asymmetric global 
coordinate rescaling such that the afiine parameter along 7 is preserved. That is, the 
Penrose limit is the metric 

ds^ = lim \-^dsl , (3.23) 
where ds\ is the metric (|2.4|) with the coordinate rescaling 

{u, V, y^) ^ (m, a V, XY") . (3.24) 



This leaves 



dr = - 2dudV + Cab{u)dY''dY'' , (3.25) 



where Cab{u) = Cab{u, 0,0) are the transverse metric components restricted to the 
geodesic 7 with V = Y"" = 0. Since, as we have just seen, these components control 
geodesic deviation and the optical tensors, it follows that the truncation of the orig- 
inal metric implied by the Penrose limit encodes precisely the information we need 
to take account of geodesic fluctuations in the context of QFT. 



The Penrose limit (|3.25|) is the metric for a plane wave in Rosen coordinates. This 



is a powerful result: it allows us with no loss of generality to consider only plane wave 
backgrounds in our analysis of photon propagation in curved spacetime. A familiar 
alternative description of the plane wave metric is in terms of Brinkmann coordinates 
{u,v,y^). Around the preferred geodesic 7 these are null Fermi coordinates [17], 
which explains their importance in our analysis of vacuum polarization. Applying 
the Penrose rescaling ( |3.23D ,( PT^^ with respect to these coordinates, the general 



metric (|3.15| ) reduces to 



ds = — 2dudv — Riuju{u)y''y^ du + S^dy^dy^ , (3.26) 
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which is the well-known form for a plane wave in Brinkmann coordinates. Note 
that this depends only on the curvature components Riuju which occur in the Jacobi 
equation. We therefore see that in both coordinate descriptions the Penrose limit 
captures the essential physics of geodesic deviation. 

The relation between the Rosen coordinates {u, V, Y"") and Brinkmann coordi- 
nates (m, V, y^) is best expressed in terms of a zweibein E\{u) which relates the 
transverse coordinates and ensures that the transverse space is flat in Brinkmann 
coordinates, i.e. 

CM = E\{u)5,,E\{u) (3.27) 

The transformation of the null coordinates is found by evaluating the geodesic equa- 
tion in Brinkmann coordinates and using the fact that in Rosen these are V = const., 
= const. This gives 

1 . . (3.28) 
V = V- -^ijyy , 

where Ei"'{u) is the inverse zweibein and we define^ 

d 

du 

Inverting, 



= Ej^—E'a . (3.29) 



y' = E\Y'' 



1„ ...... (3-30) 



with Qab = E"^ aVLijE\.^^ Now, using the definition ( |3.27|) to transform the Rosen 



form ( p.25| ) of the plane wave metric to Brinkmann form (|3.26| ), we find 



a 

hij R^njy^ Ei —j—^Ejd , (3.31) 



^fl'ij defined as Vifcj is clearly symmetric for a gradient field k^^. To prove symmetry directly 



from ( 3.29| ), note that the zweibein is parallel transported along the geodesic 7, i.e. WuE^a = 0. 



This condition implies 



rpi rpic rpk rp 

a — ^ c-t^ka ■ 



Contracting appropriately with a further zweibein then shows either Vlij = Qji or Qat = ^ba- 

^"^While flab are the only non- vanishing components of Vpfcjy in Rosen coordinates, in Brinkmann 
we find 

nf,^= \ 
-n^y n 



where on the r.h.s. fl just denotes the transverse matrix flij. Notice, using (3.34), that this satisfies 



/c^il^i, = 0; see footnote 3. 
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consistent with the exphcit result (|3]^) for the curvature. Since this also shows that 



hij 



j, we confirm that fijj defined here is the same as the previous 



du 



definition I],;, = Vj/cj. This is also immediately apparent from ( p.34| ) below. 



The eikonal ansatz for the electromagnetic field is very simple in the plane wave 
gravitational metric. In Rosen coordinates, the ansatz (|2.1|) is realized with the tree- 



level eikonal phase G = —uV, which gives = (1,0,0), while solving ( p.5|) for the 
amphtude and polarization we find 

A oc det{E\y^ , (3.32) 

^{i)a = . (3.33) 

the latter being the solution of dui{i)a ~ ^a'^{i)b = 0. In Brinkmann coordinates, 
where O = uj{v + ^yfly), we have the same result for the amplitude while, 

fc^ = {l,-yC^n + Q^)y,Q,ky') , (3.34) 
%) = {-n,ky\0,6,,) . (3.35) 



3.3 The Van VIeck-Morette matrix 



A key ingredient in the quantum field theory analysis of propagators and vacuum 
polarization in curved spacetime is the Van Vleck-Morette determinant. This is 
because it captures the geometry of geodesic deviation. It also has the important 
property of becoming singular at conjugate points along a null geodesic which, as we 
shall see, is the key to understanding the all-important analytic structure of Green 
functions. 



The Van Vleck-Morette matrix is defined from the geodesic interval 

1 /•! 

a{x,x') = -- dTg^^{x)x''x'' , (3.36) 



2^0 

where x^{t) is the null geodesic joining x = x{0) and x' = x{l), and is 

A (x x') - ^'^^"^'"^^^ (3 37) 



The VVM matrix is particularly simple in the plane wave background, primarily 
due to the isometry of the metric with Killing vector dy- First, consider the geodesic 
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interval in Rosen coordinates. We have, 



^(u-u') I du" {2V -CnhY^'Y^ 



{u~u'){V -V) - -{u-u') 



(3.38) 



where we have used the geodesic equation Y'^ — 2Q'^hY^ = or ■^(^CabY'') = 0. 
Solving this gives Y"' = C^^iu)^}, with constant (note that the original congruence 
with Y"" = constant is just the special case C,b = 0) and so 



{Y-Y' 
Substituting back, we find 



du" c^u") . 



(3.39) 



where we define 



{u - u'){V - V) - -Aab{u,u'){Y - Y'nY - Y 



(3.40) 



Aab(u, u') 



[U — U 



-\ -1 



du C 



(3.41) 



Aab{u,u') are the transverse components of the Van Vleck-Morette matrix, the de- 
terminant itself being simply 



- detAab{u,u) 
dependent only on the 2 dimensional transverse subspace. 



(3.42) 



The corresponding result in Brinkmann coordinates is similar. Here, we have 

' " \n - u') r du" (2v + h.^y - y%,y^) 

^ (3.43) 



a[x, X 



= (m - u){v -v')--{u- u) [y'SijV^jl, ■ 

In the second line, we have again used the geodesic equation + h'^jy^ = 0. Since 
this is just the Jacobi equation ( p.l7|) , we can use the expansion (|3.18|) to rewrite the 
geodesic interval in terms of the A and B functions as follows: 



a 



{x,x') = {u — u){v — v') — -{u — u) y{u)[A {u,u) ~ A ^(u ,u))y{u') 
+ y{u)A-\u',u)B{u',u)y{u) - y{u')A-\u,u')B{u,u')y{u') 



(3.44) 
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Finally, using the property A^(u,u') = ~A{u',u), we find the transverse VVM 
matrix components in Brinkmann coordinates: 

Ai,{u,u') = -{u-u')Ajl{u,u') . (3.45) 



Both these forms ( p.41| ) and ( p. 451) for the transverse VVM matrix components 



depend on quantities governing geodesic deviation. To see their equivalence, note 
that since both A and the zweibein E satisfy the Jacobi equation, their Wronskian 
is a constant matrix K, i.e. 

E^A - e'^A = K . (3.46) 

Integrating this equation, and imposing the boundary conditions defining A{u,u'), 
we find 

A,iu,u') = EM / du" C^'iu") Ejhiu') , (3.47) 

Ju' 

verifying the transpose/antisymmetry property of A{u, u') used above. This confirms 
the expected relation between the transverse VVM matrix components in Rosen and 
Brinkmann, i.e. 

Aab{u,u') = E\{u)A,,{u,u')Eh{u') . (3.48) 

This allows us to find the Rosen VVM matrix directly by solving the Jacobi equation 
with appropriate boundary conditions for A{u,u'), rather than using the integral 
form (|3.41| ). As noted at the end of section 3.1, this turns out to be the most 



efficient method of evaluation in the QFT examples that follow. 



4. Refractive Index for Scalar QED 



With these geometric preliminaries complete, we now return to quantum field the- 
ory and the calculation of the refractive index for photons propagating in curved 
spacetime. In this section, we review briefly the results for scalar QED derived in 
refs. [1-3]; our new results on the generalization to spinor QED are given in sect. 5. 



4.1 Vacuum polarization and the refractive index 

To find the refractive index at 0{a), we need to evaluate the r.h.s. of ( p.7| ) with 
A'^{x') approximated by the eikonal ansatz (O) with Q{x') = —ujV. In scalar 
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QED, the one- loop vacuum polarization tensor is given by the Feynman diagrams 
shown in Fig. |l|, viz: 



U];l°°P{x,x') = 4nagf,J^^\x - x')G{x,x') 



d^G{x, x')&^G{x, x') — G{x, x')d^dlG{x, x') 



(4.1) 




Figure 1: The two Feynman diagrams that contribute to the vacuum polarization to order a. 



The scalar propagator G(x, x') in a general curved spacetime can be expressed 
in the following form, using standard heat-kernel or proper-time methods: 



G{x,x') = ^""^^^'^^"'"Z / 7^^e-^"^'^+2#'^(^'^')^(x,x'|T). (4.2) 



^detA^^(x,a;') 1"^ dT • , i 
{g{x)g{x')Y" Jo Wry 

This arises from a representation of the propagator G{x, x') in terms of a path inte- 
gral over fluctuations around the classical path from x to x'. The VVM determinant 
accounts for the Gaussian terms, while the factor Q{x,x'\T) encodes higher correc- 
tions. It can be expanded in the form fl{x,x'\T) = 1 + J2'^=i^n(yX,x')T"' where 
the coefficients functions of the curvature and are well-known at low order. 

Clearly, this translates into an expansion in powers of R/m? in the propagator. 

Now focus on the second diagram in Fig. |l]. Inserting the eikonal ansatz in ( |2.7| ) 
and isolating the exponent term, we find 

v^rfVni-rp(x,x')A-(x') = 



(4.3) 



where we have introduced a change of variable Ti = and T2 = T(l — ^) on the 
proper time parameters for the two propagators. Now, in the limit R <^ m^, we 
can evaluate the integral over x' using a stationary phase approximation, found by 
extremizing this exponent with respect to x', i.e. 

^ d'^a{x,x') -ud'^V = . (4.4) 



2Te(l-0 



- 19 - 



Since d''^a{x, x') is the tangent vector at x' of the geodesic passing through x' and x, 
the stationary phase solution corresponds to the null geodesic with tangent vector 
QiiJ.yi — kf^i^x'), i.e. the original photon trajectory. Assuming that the background 
spacetime has an isometry we can write simply 

dv'<j{x, x') = —{u — u') , (4.5) 

and so determine u' from 

u-u' = 2uTi{l - i) . (4.6) 



This is a crucial step in our analysis. The use of the stationary phase method at 
this point allows us to go beyond the expansion in powers of uj used in the all-orders 
derivative expansion of the effective action in refs. [13, 14]. Capturing the non- 
perturbative dependence of the refractive index on lo is essential in determining its 
high-frequency behaviour, which in turn is necessary for understanding how causality 
is realized. 



The next critical observation is that the corrections to the vacuum polarization 
n^"3°°P(a;, x') at leading order in the weak-curvature R/m? expansion then arise from 
the Gaussian fluctuations around the classical null geodesic joining x and x' . It 
follows that the relevant aspects of the geometry of the background spacetime are 
simply those encoded in its Penrose limit. Moreover, the asymmetric scaling ( p.24| ) 



also ensures that the stationary phase solution (|4.6| ) remains invariant in this limit. 
So to leading order in R/m?, we only need evaluate n^"j]°°P(x, x') in the Penrose plane 
wave corresponding to the given curved spacetime and classical photon trajectory. 

From now on, therefore, we work exclusively in the Penrose limit. In this case, 
we can adopt Rosen coordinates and use the simple expression ( p.40|) for the geodesic 
interval, as well as the solutions ( p.32| ) and ( p.33| ) for the amplitude and polarization. 
From ( p.l2| ), the refractive index is 



1 - ^ / dx'^) ^ i^^ix) ni-^P(x, x') q^^ix') e^-(^-^') . (4.7) 



For convenience, we choose the origin of coordinates so that x = (m, 0,0,0). The 
integral over V is trivial and, using ( |3.40| ), leads to a delta function constraint 
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which saturates the u' integral. This automatically enforces the condition (^4.6| ), since 
the stationary phase solution becomes exact for the plane wave background. 

Since e^^-^ = (0, 0, Ei"-) is only non-vanishing in the Y"' directions, the contractions 
in ( [4.7| ) pick out just the transverse 11^^^°°'' components. Since the only dependence 
on y and Y"^ occurs in the exponents, the derivatives in ( |4.1| ) just act on these 
factors and the integrals over Y'"- are simply Gaussian. These are readily evaluated 
in terms of the VVM determinant: 



dY'- dY'^ 2 ^detA,,(n,«') ' 

(4.9) 

To complete the calculation, we now substitute the Rosen solutions for the po- 
larizations and amplitudes, and simplify the various factors involving the metric 
determinant and zweibeins. Then, the contribution of the first Feynman diagram in 
Fig. |l| is added, which removes the UV divergence at T = 0. Finally, to express the 
refractive index in its most useful form, we use ( 3.48| ) to convert the VVM factors 



from Rosen back to Brinkmann. This gives the remarkably elegant result: 

(4.10) 





X / di 





A(m, u') v^detA(M,M') 



=u-2u)T^{\~e,) 



where A is the WM matrix in Brinkmann coordinates, i.e. with elements Ajj. Of 
course, implicitly n{u; uj) is evaluated at a point x{u) on the geodesic 7. 

This is our final result for the refractive index and gives the full frequency depen- 
dence of the phase velocity of photons travelling through an arbitrary background 
spacetime. The key insight, that to this order the quantum effects on photon prop- 
agation are entirely determined by the geometry of geodesic fluctuations around the 
classical null trajectory and are therefore encoded in the plane-wave Penrose limit of 
the original spacetime, explains why the final result depends so simply on the VVM 
matrix only. 



-'^^Also note that in a plane wave background, the form (4.2) of the propagator is WKB exact 
with Vl{x,x'\T) = 1, so in this special case the leading-order results are actually exact for all R, m 
and (jj. For a general spacetime, our analysis applies in the limits R ^ and w ^ VR, the latter 
imposed by the eikonal approximation. 

This is a significant simplification over performing the whole calculation explicitly in Brinkmann 



12 



coordinates, where the geodesic interval is (3.44) and the relation 



B{u,u') + A{u,u')n{u') = E{u)E-'^^(u') 



between A and B, equivalent to (3.20), has to be used 
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An immediate consistency check is to recover the low-frequency hmit of the re- 
fractive index directly from ( [4 .101) and compare with the results of section 2 obtained 
from the effective action. The low-frequency behaviour is found by expanding the 
WM matrices in powers of t = 2ci;^(l — ^)T, since this is equivalent to expanding 
the refractive index itself in powers of uj'^R/m^. Substituting this expansion, 



Aij{u,u 



into (|;Tg), we find 

, scalar 



n. 



{u; uj) = 6ij + 



t) 



a 



Sij + -Rmujiu)t + Oif 



SGOvrm- 



R / 



'u^R 



(4.11) 



(4.12) 



in Brinkmann coordinates, in agreement with ( 2.17|) derived from the effective action. 



4.2 Analyticity 

The analyticity properties of the refractive index are essential in understanding how 
causality is realized for QED in curved spacetime. We summarize here (see ref. [1] 
for an extensive discussion) how the singularities in the VVM determinant induced 
by the existence of conjugate points give rise to a novel analytic structure for the 
refractive index in the complex u plane. This makes 5'-matrix theory and dispersion 
relations very different in curved spacetime from the familiar flat spacetime axioms 
and theorems. 



It is convenient first to rewrite ( 4.10| ) in the form 



where 



ro°-i^ rft r / 

J^{u;z)= / 1 - A{u,u ~ t) J detA{u,u - t) 



(4.14) 



Notice that we have to introduce a prescription, as indicated, for dealing with the 
branch-point singularities of the integrand that arise whenever x{u) and x{u — t) are 
conjugate points and det A(-u, u — t) diverges. This prescription actually follows from 
a careful treatment of the VVM determinant factor in the propagator {e.g. see the 
book [18]). Where we have written det A^i,{u,u — t), it should be interpreted as 

e-^'^'/^ldet A^^(M,M-t)| , (4.15) 

where u is the Maslov-Morse Index which counts the number of conjugate points 
(or more properly the number of times an eigenvalue of A^y{x,x') diverges) on the 
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geodesic joining x{u) and x{u — t). Another way to interpret the phase is that 
it provides the prescription for the analytic continuation of A^i^{u,u — t) into the 
complex t plane. In particular, it requires the t contour to avoid the singularity at 
a conjugate point by veering into the lower-half t plane as indicated above. This 
choice is consistent with the requirement that the flat-space limit is smooth. Note 
that since the integral is over Ret > 0, it receives support only from that part of the 
null geodesic to the past of i.e. from A{u,u — t) with Ret > 0, as expected in 
a causal theory. 

The most important observation is that ^{u; z) is by construction guaranteed 
to be analytic in the lower-half z plane (including the positive real axis). Since z 
and u are inversely related this means that n{u;u!) is analytic in the upper-half of 
the u complex plane (including the positive real axis). This establishes the first 
requirement for micro-causality. The second requirement rests on the behaviour of 
the large |ti;| limit of the refractive index. Large \uj\, corresponds to z — ^ 0, in whose 
limit 



The question here is whether the t integral defining ^{u; 0) is convergent. If the 
spacetime is asymptotically flat in the far past, then A.{u,u — t) asymptotes to a 
constant and the integral is indeed convergent. However, this is sufficient but not 
necessary as we shall see with particular examples. 

Assuming that n(u; uj) is analytic for large uj in the upper-half plane then we 
can write the Kramers- Kronig relation in the form ( |1.2|) , where the contour avoids 
any non-analyticity on the real axis but lying just above the real axis. What we shall 
find is the condition n{—uj + ie) = n{uj + ie)* is not satisfied in any of the curved 
space cases. This should not be surprising since the background geometry breaks 
Poincare invariance. 



5. Refractive Index for Spinor QED 

In our previous work [1-3] we have focused on scalar QED in order to investigate 
the novel physics of photon propagation in curved spacetime in a relatively simple 
context. Here, we develop the extra formalism needed to study the realistic case of 
spinor QED and derive an explicit expression for the full frequency dependence of 
the refractive index in this theory. 




(4.16) 
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5.1 Spinor propagator in a plane wave spacetime 



In order to define spinors in curved spacetime, we first introduce a local pseudo- 
orthonormal frame at each point in spacetime by means of a vierbein: 



where we preserve the null structure in the local frame by taking 

/ -1 o\ 
■10 
10 

ly 

For a plane wave spacetime, in Rosen coordinates, the vierbeins are explicitly 

1 

0E\ 



(5.1) 



V 



(5.2) 



(5.3) 



In the local frame we denote the indices as A = (+, Notice that the indices 

i,j — 1,2 are common to the local frame and the Brinkmann coordinates in the 
transverse space. The 7 matrices are defined in the local frame and satisfy 



{7^,7^} = -2v 



AB 



(5.4) 



In particular, along the null directions the gamma matrices are nilpotent, 7+7+ = 
7~7~ = 0. This is a crucial simplification, as we see below. 

The next step is to define the spin connection. From the metric condition 
^n^-^B = 0, this is: 

(^fiAB = GAud^.e^B + eAu^" ^ipe'' B , (5.5) 
and has non- vanishing components 



^oi+ — — ^a+i — — Eia{u) . 



(5.6) 



One can verify explicitly that the connection is torsion free: dcA + i-oab A = 0, or 
in components, 

d[^,eAu] + i^[f,ABe^u] = . (5.7) 
The covariant derivative on spinors is defined as 



where a^^ — 7 [7"^, 7^]. In Rosen components, this gives 



(5.8) 



(5.9) 
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so the Dirac operator is 

y = 7^e^^V^ = j+du + rdv - fE,%u)da + \E,\u)E^,{u)^y^+ . (5.10) 

The spinor parallel transporter is a bi-spinor that parallel transports a spinor 
along a given path, in our case the null geodesic 7 joining x and x' . It satisfies the 
two conditions 

d^a{x,x')V^\J{x,x') = Q , U(x,x)=I, (5.11) 

and has an explicit representation in terms of the spin connection by the path ordered 
expression 

U(x,x') = Pexp-- / uj^Asa^'^ dx^ , (5.12) 
2 Ju' 

where the integral is taken along 7. The key observation is that since the spin 
connection involves the nilpotent 7"'', we can expand the exponential and only terms 
up to linear order in the spin connection contribute (in particular the issue of path 
ordering is moot): 

l-Y-{u) 

U(x, x')=l+ IY7+ / Eia{u) dY'' . (5.13) 
To evaluate this, write 

/ EUu)dY''= duEia{u)Y\u) . (5.14) 

JY"^ Ju' 

and, as in section 3.3, use the geodesic equation = C°*^f,. Since C = E and 
f2 = EE~^ are both symmetric, we can write 



j duEC = JduEE-\E-^y= J du {E-^y E^ (E-^) 

= - j du{E'y = -{E-y . 



(5.15) 



Hence 



I EMdY'^ = -{E,'^{u)-ErW)% 

= - {E,\u) - Eau'))^^^^^^{Y - Y'Y 



(5.16) 



u — u 



using ( p.39| ) and ( p.41|) . This gives the following explicit form for the spinor parallel 



transporter in terms of the VVM matrix: 

U(x,x') = I - \^^^^{E,\u) - E,\u'))^^!^!^^{Y - Y'f . (5.17) 



u — w 
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Finally, we need the Feynman propagator S{x, x') for the spinor electron. This 
can be written in terms of the propagator G{x, x') for a massive scalar field and the 
spinor parallel transporter \]{x, x') as follows: 

S{x, x') = - m)G{x, x')l]{x, x') . (5.18) 



5.2 Vacuum polarization in spinor QED 



With this form for the spinor propagator in a plane wave spacetime, we can now cal- 
culate the vacuum polarization and refractive index following the method of section 
4. The vacuum polarization is given by the usual one-loop Feynman diagram (the 
second in Fig. p, viz. 

n;-j°°P(x,x') = e'Trlf,^ix,x'h,Six\x)] . (5.19) 



Inserting the vacuum polarization for spinor QED into (|2.12|) for the one-loop 
contribution to the refractive index, we find, analogously to (|4.7|) in the scalar case, 

n,,{x;uj) = 6,, + % f rfV v/^^TT Tr [7.^(2;, x')7,S(a;', x)] e^-(^-^') . (5.20) 

Again, the calculation is best performed in Rosen coordinates although the end result 
is most naturally expressed in terms of tensors with Brinkmann indices. We again 
take X = {u, 0, 0, 0) for convenience. 



Now notice that the terms which are linear in m involve an odd number of 
gamma matrices and so do not contribute to the trace. There are two remaining 
contributions to the trace of the form 



Tr 



-fif (G(x, X )U(x, x')) -fjf {G{x', x)U{x', x)) 



+ m Tr 



7iG(a;, x')U(x, x')'yjG{x' , x)lJ{x' , x) 



(5.21) 



The explicit expression for the spinor parallel transporter is given in ( |5.17| ). The 
integrals are Gaussian and easily evaluated while the V integral is trivial and simply 
produces a delta function constraint which enforces the same condition ([4.6|) as in 
the scalar case. 



Noting that 



yU(x, xOy' U(a;', x) = 



(5.22) 
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due to the fact that 7+7+ = 0, there are three types of contribution in ( 5.21|) which, 
suppressing the calculational details which are not in themselves enlightening, yield: 



a X \/g[x') e . , , Tr 



A{x 

= 647rT2^2(l - OVdetA(u,u' 
along with 

A{x' 



'yiG{x, x')V{x, x')'jjG{x' , x)V{x' , x) 



(5.23) 



d^x' \/g{x') ' 



A{x] 



Tr 



Ii1/G{x, x')U(x, x')7jy'G(x', x)U(x', x) 



mnTiil - OVdetA(M,M'){ ( - 8^(1 - 



- (1 - 2^(1 - - w')Tr (5,' A(n, u')^-\u, u')))5,, + 4^(1 - 0A.,(«, W 



(5.24) 



and 



a;/)e-*^^'4^Tr 



A{x) 



-fifG{x, x')U(x, a;')7j- G(x', x)f'l]{x', x) 



+ 7j G(x, x'Jf\]{x, x')'^j^'G{x', x)U(x', x) 



(5.25) 



■U, M' 



= 16mT^(l - i)5ij^deiA{i 

X Tr (^1 - A(m, m') + (m - u')du'^{u, u')^-^{u, m')) • 
In all these expression u' is constrained via 



Putting everything together, the final result for the refractive index is 
a r dT 



27rcu2 ./n r2 



■ le 



X 



^ v/detA(M,M') (^4iTm2 + ^^^^-^Tr (A(m, m') - 1 
2{u - u)Tt {du'A{u,u)A-\u,u')) + 8)1 - 4A{u,u' 



J n'=M-2ajT5(l-C) 



(5.26) 



or equivalently, 



X 



a 



dT 



It 



a2T 



A/detA(M,M') 



e(i-0 



Tr(A(u,u') - 1) 



- 4(m - u')Tr (9„' A(u, m')^" K u')) + 4 1 - 4A(m, u' 



J u'=M-2a;T5(l-0 

(5.27) 
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In the second expression, we have taken the first term and integrated by parts in T 
and ignored the singular boundary term which is curvature independent. The result 
is then manifestly UV finite, i.e. there is no singularity of the integrand at small T. 



Finally, changing variables to t = 2u;^(l — and writing the refractive index 
in the form of section 4.2, we have 



"••■'■■"(";-) = i-£;/<««i-0^''""('';^:;|;^) . (5.28 



where 



- ze"*^Vdet A(«, u-t) [(^^^-^Tr (A(«, u-t)-l) 
+ 4iTr(9fA(M,M-t)A-^(M,M-t)) +4)1 -4A(M,M-t) 



(5.29) 



The generic analyticity properties are exactly as described for scalar QED. Also 
as in the scalar case, we can verify the leading low-frequency term in the expansion 
oin{u]uj). Substituting the expansion ( [4.13|) of the VVM matrix into ( p.29|) , we find 



\u-u;) = 5,,-Y^£^(m,i?„„-4i?_,) + §iO[^) , (5.30) 



spinor / ^ <- ^ / ^ ^ ^ . ^ \ R / ^ 

' — — ' -11 ') \ Lfc u, u,t,u, 1 / • ') ~ \ 

in agreement with the original Drummond-Hathrell formula (p.l8|) 



6. Homogeneous plane waves 



The first class of examples we shall study in detail are the homogeneous plane waves. 
These are already sufficient to display many of the interesting consequences of curved 
spacetime geometry on the analytic and causal structure of the refractive index. 
However, as we shall see, they are not just simplified toy models but actually arise 
as the Penrose limits of physically interesting spacetimes. 

Homogeneous plane waves are characterized by an enhanced degree of symmetry. 
Any plane wave admits a Heisenberg algebra as its spacetime isometry (though this 
is not necessarily a symmetry of the original metric for which the plane wave is a 
Penrose limit.) However, there are two simple cases (see [12] for a more complete 
classification) where this isometry group is extended - the symmetric plane waves, 
with profile function hij constant, and the singular homogeneous plane waves, where 
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hij = Cij/u^. The first case is evidently invariant under it-translations; the second 
under scale transformations ti— >Att, v ^ \~^v. 

An important insight in the study of Penrose limits is the idea of hereditary 
properties, i. e. some feature of the original metric which is preserved in the Penrose 
limit [12]. Amongst such hereditary properties are (i) Ricci flat (which implies 
tr/ijj =0), (ii) conformally, or Weyl, flat {hij — h{u)5ij, so the Penrose limit is 
cylindrically symmetric), (iii) locally symmetric, i.e. the covariant derivative of the 
Riemann tensor vanishes. Another useful property is that an Einstein metric has a 
Penrose limit which is Ricci flat. 

Consider first maximally symmetric spacetimes, including de Sitter (with isome- 
try S0(l,4)) and anti de Sitter (SO(2,3)). The symmetry implies that the Riemann 
tensor has the form R^i^xp — | [gixxgvp — gnp9ux) , which implies this is both an Ein- 
stein space and conformally fiat. Since this requires the Penrose limit (which is 
unique because the maximal symmetry implies all geodesies are equivalent) to be 
both Weyl and Ricci fiat, hij must be proportional to 6ij and traceless, and therefore 
vanish. That is, the Penrose limit of a maximally symmetric spacetime is fiat.^^ 
Physically, this means that the refractive index receives no quantum corrections for 
photon propagation in de Sitter or anti de Sitter space. 

Next, consider locally symmetric spacetimes, i.e. with V aR^vXp — 0. Since 
this is an inherited property, the Penrose limit has a profile function hij which is 
independent of u. This defines a symmetric plane wave. With no loss of generality, 
we can diagonalize hij and write the metric in Brinkmann coordinates as 

ds^ = -2dudv-hijyydu'^ + dy'dy\ (6.1) 

where hij — diag((Ti,(j|) with the constant. The corresponding curvatures are 
then Ruu — erf + and Ruiui — ■ There are two generic cases to consider: (i) ctj 
both real (which includes the conformally fiat case ai — and (ii) cxi real and 
(72 = i\o-2\ pure imaginary (including the Ricci fiat case ai = \(J2\)- The third case, 
with both a I and (72 imaginary is disallowed by the null energy condition, which 
requires Ruu > 0. 

^^We can derive these properties very directly in the Newman-Penrose formahsm. (See refs. [5,6, 
9, 19] for an extensive discussion of the NP formalism applied to superluminal photon propagation 
in curved spacetime.) As explained further in section 7, the NP tetrad provides a realization of the 
local null Fermi coordinate basis, with the Penrose limit profile function being determined by Ruiuj, 
that is by R^^xpi^^m" i-'^m'' and Rp_^xpi'^m''£^fh'' , or equivalently Cimem = — *o and i?« = — 2<i?oo- 
For a maximally symmetric space with Rf^^xp = l{gp.xg„p - g^pQi^x), Rtmtm is proportional to 
(^rn? — £.m£.m (similarly for Rem£rn), which vanishes by the defining properties of the NP tetrad. 
Both $00 and \E'o vanish and the Penrose limit is flat. This formalism also makes it clear why 
the Penrose limit of an Einstein space is Ricci flat: if i?^^ = A^^^^ in the original space, then 
$00 = -iRuu = -IR^J^P" ~ ^2 = 0. 
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Symmetric plane waves arise as the Pemose limit for critical photon orbits in 
Schwarzschild spacetime, as shown in section 7, and also in higher dimensions as the 
Penrose limits corresponding to a class of null geodesies in the AdSs x spaces 
of interest in the AdS/CFT correspondence. The Penrose limit for geodesies lying 
entirely within the AdS or subspaces are flat, since these subspaces are maxi- 
mally symmetric. However, a null geodesic in AdS with a non-vanishing angular 
momentum in the subspace has a non-trivial Penrose limit. Since the original 
space is locally symmetric, and this is a hereditary property, the Penrose limit for 
these geodesies is a symmetric plane wave and in fact [12] can be shown to have 
hij ~ diag((j^, . . . , erf , cr|, . . . cr|), with four factors of af and p — 1 of cr|. (Here, 
ri = crj~^ and ^2 = cTg"^ are the curvature radii of the AdS and subspaces respec- 
tively.) 

The second class of homogeneous plane waves considered here are the singular 
homogeneous plane waves, with Brinkmann metric 

ds^ = — 2dudv — Cijy'^y^ — — + Sijdy^dy^ , (6.2) 
where, again, with no loss of generality we can take Cj,- = ( ^ ^ ] to be diagonal. 



C2, 

Here, we find that the singularity structure of the VVM matrix, and consequently 
the analytic properties of the refractive index, turn out to be dependent on the 
actual numerical values of the constants Cj. These metrics arise as the Penrose limits 
associated with spacetime singularities, with the values of q being related to the 
Szekeres-Iyers [20,21] classification of power-law singularities [12]. In particular, the 
singular homogeneous plane waves arise as Penrose limits both of cosmological FRW 
spacetimes and of the near-singularity limits of black hole spacetimes. 



6.1 Symmetric plane waves 

Now consider the symmetric plane waves, with Brinkmann metric ( |6.1| ). The Jacobi 
equation ( ^.IT] ) for the connecting vector ?/* is simply 

-^y' + afy' = , (6.3) 

with general solution y^{u) = a cos{aiU + b). The VVM matrix is then found from the 
matrix A{u,u') which solves the Jacobi equation ( |3.19| ) subject to the appropriate 
boundary conditions. This selects the solution 

A,,{u, v!) = 6ij '-^ . (6.4) 
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so we find 



, . ,^ . (Tiju - U') 

Aij{u, u ) = dij- - 

smai[u — u') 



(6.5) 



Many of tlie key features of the analytic structure of the refractive index are 
already evident in the simple case of a conformally fiat plane wave, cxi =02 = o . 
Recall that the refractive index is given in terms of the auxiliary function JF(m; z) 
by ( 4.13|) . For scalar QED, inserting the formula ( |6.5|) for the VVM matrix into the 
definition ([4.14| ), we find J-{u\ z) = J-'{z)l with 



scalar / 



dt 



-izt 



at 



sin at 



(6.6) 



while for spinor QED, ( ^.29|) gives 



spinor ; 



Z) 



dt 



le 



-izt 



at 



sin crt 



2 cos crt + 1 

2{3-at ) + 

sm at 



^(1 -0 Vsinat 



at 



(6.7) 

The integrands here have a series of poles (a special property of the conformally fiat 
case; in general these will be branch points) on the positive real axis at t = nvr/s, n = 
1,2, .. . corresponding to the singularities in the VVM matrix related to conjugate 
points on the classical null geodesic [1-3]. 

The refractive index itself is then evaluated by performing the integral over 
t, with the appropriate contour rotation described in section 4. In this case, the 
refractive index is real, Imn(ci;) = 0, since the contour can be rotated t —>■ —it to 
lie along the negative imaginary axis in which case the integrand is manifestly real. 
In both cases the t integral can be performed explicitly to give, for scalar QED (a 
result from [1])^^ 



^scalar/ N 



:log(^)+zV^(l + ^)-a 



(6.^ 



while for spinor QED 



jTspinor^^^) = 8( - (T - 3z + 3alog(27r) - 6(Tlogr(| + ^) 
+ ^^(l + ^) + 2^^(i + ^)) 

4 .(a + 2;-alog(27r) + 2alogr(i + ^)-#(l + ^)) . 



(6.9) 



+ 



e(i-o 



In both cases J-'{z) is analytic in the lower half plane (including the positive real 
axis); in fact the only points of non-analyticity are ai z = 0, —a, —2a, . . ., being 
either branch points or poles, or combinations thereof. 



14 



Here, V'(^) = dzT(z) is the di-gamma function. 
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Figure 2: The refractive index n(Lu) — 1 of a conformally flat symmetric plane wave, in units of 
acr^ /(27rTO^), plotted as a function of logwcr/m^ for scalar QED (left) and spinor QED (right). 

The remaining ^ integrals for both scalar and spinor QED can be performed 
numerically and the results for the refractive index are plotted in Fig. This shows 
clearly how the occurrence of a superluminal low-frequency phase velocity in spinor 
QED is compatible with causality, n{oo) 1. Also note that since lmn{uj) = 0, the 
Kramers-Kronig dispersion relation in its conventional form ( |1.1| ) cannot hold since 
it is clear that n{uj + ie) ^ n{~uj + ie)* for uj G R^. However, the more general 
dispersion relation ( |1.2|) does hold as demonstrated in full detail in ref. [1]. 

Another interesting example is the Ricci flat symmetric plane wave, ai = a and 
(72 = icr. In this case, the VVM matrix has components 



with corresponding expressions for ^{z) for scalar and spinor QED. The t integrands 
therefore have series of branch points on both the real and imaginary axes. As a 
result, when the rotation of the t contour, t —>■ —it is made so that it lies along the 
negative imaginary axis, we encounter the new branch points, which contribute a 
non- vanishing imaginary part to the refractive index. (See ref. [1] for full details and 
discussion of the analytic structure of ^{z).) 

The numerical results for the refractive index are plotted in Fig. ^. Consistent 
with the low-frequency results ( p. 171) , (|2.18| ), the two polarizations give opposite sign 
corrections to the refractive index (gravitational birefringence), in contrast to the 
conformally fiat case. The imaginary parts are both positive for spinor QED, as ex- 
pected from the fiat spacetime identification of imaginary parts of forward scattering 
amplitudes with cross sections (optical theorem), but are both negative for scalar 
QED. This occurrence of negative imaginary parts of the refractive index, which is 
repeated in many of the examples which follow, remains to be fully understood. 





(6.10) 



- 32 - 



0.004 










0.002 

















--" 6 


i 10 


0.002 




/ 














0.004 




7' 










/ 






0.006 




/ 








Figure 3: The refractive index n{uj) — 1 of a Ricci flat symmetric plane wave, in units of 
a(T^/(27rm^), plotted as a function of log a;cr^/(27rTO^): continuous (real part, polarization i — 1); 
big dashes (imaginary part, i = 1); small dashes (real part, polarization i = 2) and dots (imaginary 
part, i — 2), for scalar QED (left) and spinor QED (right). 

6.2 Singular homogeneous plane waves 



We now come to the singular homogeneous plane waves, with Brinkmann metric 
( |6^ . Here, the Jacobi equation is 



y + = 



du 



(6.11) 



and it is convenient to define at = \/Y^^^Ci so that the general solution is = 
(jy^ii+a,)/2 _|_ /}y(i-a»)/2_ select the (diagonal) matrix Aij, we impose the familiar 



boundary conditions A{u',u') = 0, -^A{u,u')\u=u' = 1 which gives 



The VVM matrix A.{u,u') therefore has components 

u — u' 



1,2. 



— U 



icii 



1,2. 



(6.12) 



(6.13) 



Notice in this case, there are no conjugate points even when is real. However, 
when Ci > |, ttj purely imaginary there are an infinite sequence of conjugate points 
when 

u = UQ^^\—- \ , n = l,2,.... (6.14) 
V fti / 



Substituting into (|4.1(]| ) for the refractive index in scalar QED determines the 
function J-{z) = J-'(^i){z)5ij as 



dt 



—izt 



1 — A'^(^(m, M — t)A22^(M, M — t) 



1/2, 



(6.15) 



for the polarization corresponding to z = 1 (similarly for i = 2), with the related, 
more complicated form (|5.28|) for spinor QED. However, the expression above only 
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makes sense when the singularity at m = hes in the future of the point u. In other 
words M < 0. This is what happens in the case of black hole as we describe in the 
next section. On the other hand, for a cosmological spacetime the singularity occurs 
in the past of the point u. In this case the expression above does not make sense 
unless the t integral is cut-off at some point uq > 0. This will discussed in more 
detail in Section |. 

In the former case, we can easily calculate the first few terms in the frequency 
expansion of the refractive index riij = riiSij-. 

n,{u- I ^^^^'^^ ^. au;(3ci + C2) ^ 

with a similar expression for the other polarization. 



7. Black Holes and their Singularities 



After studying these especially simple examples, we now turn to spacetimes arising 
as physically important solutions of the Einstein field equations and begin with black 
holes. This immediately raises the question of the effect of horizons and singularities 
on the causality and analyticity properties of photon propagation. The first step is 
to determine the Penrose limit relevant for various classical trajectories. It turns out 
that we can do this in considerable generality for black hole spacetimes. Remarkably, 
the near- singularity limits turn out to be just the singular homogeneous plane waves 
discussed above. 



7.1 Penrose limit for Petrov type D metrics 



The relation between the Penrose limit and the geometry of geodesic deviation ex- 
plored in section 3 provides a very efficient way of determining the profile function 
hij of the equivalent plane wave. The starting point is the identification of the par- 
allel transported frame Ea introduced in ( |3.12| ) in the determination of null Fermi 
coordinates. The component Eu'^ = is simply the tangent vector along the null 
geodesic 7; in fact, we then only need the explicit form for the transverse spatial 
components Ei, i = 1,2 (related to the subspace V of footnote 3). Assume now that 
7 is just one integral curve of a known vector field k^{x), or equivalently assume a 
particular embedding of 7 into some null congruence. Having identified the trans- 
verse subspace, we can immediately compute the now familiar matrix Qij = Vikj, 
first introduced through the Lie derivative (|3.2|), which defines the optical tensors 
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controlling geodesic flow. The proflle function of the Penrose plane wave metric, in 
Brinkmann coordinates, is then given by the identity hij = Qij + flik^'^j- 

In section 7.2, we carry through this construction explicitly for general planar 
orbits in the Schwarzschild metric. First, we present an alternative derivation which 
exploits the Newman-Penrose description of the whole class of Petrov type D space- 
times [11], which includes the Schwarzschild, Reissner-Nordstrom and Kerr black 
hole solutions. This makes very clear how the plane wave profile function depends 
on the conserved integrals of motion characterizing the classical trajectory 7. 

The Newman-Penrose (NP) formalism uses a null tetrad Ea^ = {L^, N'^, M^, M^), 
where the null vectors are defined so that L ■ N — —1, M ■ M — 1, with all other 
scalar products vanishing. The ten independent components of the Weyl tensor are 
represented by the five complex scalars: 



*0 = 














- C^^xpL'^N'^M^NP 


^4 = 


- C^^xpN^M'^N^Mf 



and we also encounter the scalar $00 = —\Rll from the set describing the Ricci 
tensor. The special feature of Petrov type D spacetimes is that we can choose a 
special NP basis n^, m^, fn") where is tangent to the principal null geodesies 
such that the only non- vanishing component of the Weyl tensor is \l/2- (We reserve 
the un-hatted notation for the components with respect to the principal null geodesic 
basis.) 

To apply this formalism here, we note that along a given null geodesic 7 with 
tangent vector = k^, we can identify the NP frame with that describing the null 
Fermi coordinates. (Of course, in the NP basis the two transverse spatial compo- 
nents are complex linear combinations of those in the Fermi basis.) Determining the 
Penrose limit then simply involves finding the components Rlmlm and Rlmlm of 
the Riemann tensor, or equivalently Rll and Clmlm (see below), since this gives 
the profile function of the Penrose plane wave. 

It follows immediately that for photons following a principal null geodesic, the 
Penrose hmit is fiat, since by definition ^0 = 0. This shows immediately that at 
C(i?/m^), for all frequencies, there is no quantum correction to the refractive index 
for photons following radial geodesies in Schwarzschild spacetime, with a similar re- 
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suit for the corresponding principal null directions in the Kerr or Reissner- Nordstrom 
metrics. 

To extend this to an arbitrary classical trajectory, we first expand the corre- 
sponding tangent vector in terms of the NP basis adapted to the principal null 
geodesies, i.e. 

= at + I3n^ + + 7*m^ , (7.2) 

and let 

= At + Bn^' + Cm^ + Drfi^ . (7.3) 

The null conditions imply a/? -77* = 0, AB-CD = and A(3 + Ba-CY-Dl = 0. 
The Weyl tensor for a Petrov type D metric can be written as (see ref. [11], chapt. 1, 
eq. (298), but note a typographical error) 

C^uXp = \E'2 1^'^^^] -[mmmm] + [inmrh] + [imnm]^ + c.c. , (7.4) 

where the notation [ . . . ] indicates the linear combination of the given vectors with 
the symmetries of the Weyl tensor. A straightforward calculation now shows^^ 

Clmlm = ^ Kl ^| , (7.5) 
where the complex quantity Kg is defined as 

Ks = 2%hL-eM -n- L-mM -m) 

(7.6) 

= 2%^AP-CY) . 



The significance of this result follows from a theorem of Walker and Penrose [10], 
stated and proved in ref. [11] (theorem 1, chapter 7), that Kg is conserved along the 
null geodesic with tangent L^, i.e. 

L'^V^Ks = . (7.7) 

This shows very directly how the Penrose limit is determined not just by the curvature 
but also by an integral of motion along the chosen geodesic 7. We shall see explicitly 
how this is realized in the examples below. 

^^By direct calculation we find 

Clmlm = -6((A/3 - C^f ^2 + (A/S - Djf 

However, the null conditions can be manipulated to show that {Ap — Cj*){A(3 — D7) = 0, proving 
that at least one of the coefficients above vanishes. This establishes eq. (|7.5| ). In the Schwarzschild 
example below, we see explicitly that (A/3— C7*) is non-zero and gives Kg, while {A(3—D^) vanishes 
as required. 
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The plane wave profile for the Penrose limit, allowing for a spacetime with a 
non-vanishing Ricci as well as Weyl tensor, is 



/j. . — 1 I i^LMLM + RlMLM + '^RlMLm) — i{RLMLM — RlMLm) 

2 \ -i{RLMLM - Rlmlm) - {Rlmlm - Rlmlm " '^Rlmlm 

(7 

From the definition of the Weyl tensor. 



(7.9) 



RfiuXp C^pXp -\- 2 (^9fj.xRup 9fj.pRvX 9u\Rpp ~l~ 9upRp,\) 
~ g {9p.\9up - 9fj.p9u\) 
and using the identity Cimlm = 0, we immediately find 

Rlmlm = Clmlm, Rlmlm = 2^ll ■ (7-10) 

The eigenvalues of hij are therefore 

h± = -jRLL ± \Clmlm\ 
= -l-oo ± l^ol , 

the latter form emphasizing again the importance of the NP scalars in the description 
of the Penrose limit and photon propagation [19]. 



Using ( [7.5|) , we have therefore shown that the Penrose limit for a general type 
D spacetime has a profile function h^j with eigenvalues 

h± = \Rll ± \\Ks?\^2\^ . (7.12) 
Note that, as expected, h^j is traceless when the original metric is Ricci fiat. 

In fact, this result can be simplified still further. Following Chandrasekhar [11], 
we can prove using the orthogonality and null properties of the L, A^, M, M basis 
(including the result that if A[3 — C-y* ^ then A(i — = 0; see footnote 15), that 



So we can finally write 



\Ks\^ = 4\'^2\~^aP ■ (7.13) 



^Rll ± 6\'^2\L-iL-n . (7.14) 



This shows very clearly that in the case of Type D spacetimes, the Penrose limit 
depends only on the tangent vector = itself; we do not even need to compute 
the remaining basis vectors in the set Ea- 
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The advantages of this method of determining the Penrose hmit are now clear: it 
only involves knowledge of the classical null geodesic 7 itself and does not invoke any 
particular embedding into a null congruence; it makes the relation with the conserved 
quantities characterizing the classical orbit manifest; and it allows the simplifying 
power of the NP formalism and Petrov classification to be exploited fully. 



7.2 Schwarzschild black hole 

The Schwarzschild metric is 

ds^ = -^dt^ + '^dr^ + r''(de^ + sin" ed(f)^) (7.15) 
A 

where A = — 2Mr. We consider planar classical orbits, and with no loss of 
generality choose = const. The first integrals of the null geodesic equations give 



t = E-, r = E^l-^, 9 = E-, = 0, (7.16) 

where E is the energy and b = L/E is the impact parameter, with L the conserved 
angular momentum. We can therefore identify the tangent vector as E, ^,0^ 
in the (t, r, 9, 0) coordinate basis, where F = a/1 — b'^A/r'^. 

To derive the Penrose limit, we first identify a parallel transported frame E^ 
along this null geodesic, as described above. We immediately have 

The transverse vierbeins E^ are then determined by the orthonormality conditions 
q^uEa^Eb^ = Vab, with tiab as implied in (|3.12|) , supplemented by the parallel 



transport condition k^'^ ^Ei^ = 0. The orthonormality conditions give 

1 4 2 

and 

^ (0, 0, 0, 1) , (7.19) 



r sin 6' 

showing explicitly that these algebraic conditions determine Ei^^ only up to the ad- 
dition of a piece proportional to k'^ itself. The remaining freedom, parameterized 
by the function c(r) above, is determined by the parallel transport condition, which 
gives c(r) = ^ J dr F~^. Note that we do not require an explicit expression for E^^^ 
for this construction. 



- 38 - 



The next step is to determine the matrix fij^, defined as 

% = Ei>'E/\/^K . (7.20) 

This calculation can be simplified by noticing that since k'^'S/ Jt^ = by the geodesic 
equation, and V Jiy is symmetric because k'^ is a gradient flow, we can add a multiple 
of A;^ to the vierbeins Ei^ in (|7.20| ) without affecting the result for VLij. In practice, 
this allows us to choose the function c(r) in (|7.18| ) as we like in order to make the 



subsequent calculation simple. Natural choices include c = F [6] and c = 
[12]. Exploiting this freedom, we find after a straightforward calculation that Vtij is 
diagonal in this basis with 

1 f Mb^\ ^ F b ^ , , 

fill = — 1 ^ , fi22 = - + ^cot^. 7.21 



Alternatively, using the geodesic equations (|7.16| ), we can express these in the more 
enlightening form: 

fill = -^log(rr), fi22 = :^log(rsin^) . (7.22) 

au au 

These results are clearly consistent with the trace identity [12] 

trfi = t^^^'Ea^Eb^V^K = V^P = ^9^(v/^P) . (7.23) 

The profile function hij = fijj + Qik^l'^j for the plane wave is easily computed from 
( ff:22D and we find 

1 1 d^ 

= -— ^22 = — — (rsin^), (7.24) 

rr du^ r sm du^ 

that is, 

hu = ^, = —T- . (7.25) 

As anticipated, tr hij = as the Penrose limit inherits the Ricci fiat property of the 
original Schwarzschild spacetime. 

In the Newman-Penrose method, we first note that the NP null tetrad for the 
Schwarzschild metric is [6,11]: 

r = i(r^ A, 0, 0), n^ = -^(r^ -A, 0, 0), = -^{0, 0, 1, -^) , 

A 2r^ ^^/2r smd 

(7.26) 

and the only non- vanishing component of the Weyl tensor is \&2 = —M/r^. We can 
therefore expand the tangent vector = k^ for the classical orbit in the NP basis 
as follows: ^ 

L^" = -(1 + F)r + ^(1 - F)n^ + -^(m^ + m^) . (7.27) 



2' ' A' ' V2' 



r 
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The circular polarization vector M^, which satisfies L.M = 0, M.M = 1 and is 
parallel transported, is identified with the vierbein found above: 



-I 4 2 

/ T T 1 

Expanding this in the standard NP basis we have 



(7.28) 



= -i-^(-i(l-c)(l + F)r + ^(l + c)(l-FK)+i(l + c)m^-i(l-c)^^ 

(7.29) 

This fixes the coefficients a, . . . , 7 and A, . . . ,D of eqs.( |7.2|) and dTTsf) . Once again 
we can use the fact that adding a multiple of to iVf^ affects neither the evaluation 
of Kg (since L'^ is null) or Clmlm- So to perform the calculation efficiently, we can 
make any convenient simplifying choice of c(r) in ( [7.29 ). From ( [7.6|) we then readily 
find the conserved Walker-Penrose quantity Kg-. 

Kg = ~V2{-M)-h. (7.30) 

We see directly that in this simple case the conserved quantity Kg for a given photon 
energy is determined by the angular momentum of the orbit. This confirms the ad- 
vertised link between the conservation laws governing the classical null geodesic and 



the Penrose limit. The Penrose plane wave profile itself (|7.25|) follows immediately. 



We can make a number of immediate observations based on these results (see 
also ref. [12]): 

(i) The explicit dependence of hij on the impact parameter b implies that hij = 
for classical null trajectories with vanishing angular momentum. This confirms that 
the Penrose limit is fiat, and therefore the refractive index receives no quantum 
corrections, for purely radial geodesies which of course are the principal null geodesies 
of the Schwarzschild metric. 



(ii) There is a closed photon orbit at a critical value of the radial coordinate r. 
Using ( |7.16|) we have f = E^(l — ^), fixing Tc = 3M. Then, requiring r = 
determines the critical impact parameter be = 3a/3M. For this orbit, the Penrose 
limit has 

which describes a Ricci fiat symmetric plane wave. The refractive index for this case 
has already been discussed in section 6. 
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(iii) In the near-singularity limit, the radial geodesic equation becomes 



- - ^ ^ 0(1) . (T.32, 

with solution = '^Mh'^v?. The Penrose limit therefore has 

= -|l. (7.33) 

Note that this is independent of the impact parameter 6, though remember the 
derivation is only valid for orbits with 6 7^ 0. Also, in order for an incoming photon 
to reach the singularity, the impact parameter must be sufficiently small, h < 6crit = 
3\/3M. The Penrose limit in the near-singularity region of Schwarzschild spacetime 
is therefore a Ricci flat singular homogeneous plane wave, with coefficients c^'^ = 
=F6/25. As explained in detail in [12,22,23], this specific coefficient is related to the 
power-law nature of the singularity for this class of black hole spacetimes. 

Having established the Penrose limit, we now need to determine the Van Vleck- 
Morette matrix A.{u,u'). This entails solving the Jacobi equation for A{u,u' with 
boundary conditions A{u',u') = 0, ■£^A{u,u')\u=u' = 1- Recall that the Jacobi 
equation for the geodesic deviation vector in Brinkmann coordinates is in general 



-y' + K^y^ = , (7.34) 



du 

and in the present case, hij is diagonal. Beginning with i = 1, note that the special 
form ( |7. 24| ) of the expression for hu shows immediately that a particular solution is 



y^ = fi{u) = r{u)r{u) with r{u) the solution of the geodesic equation (|7.24|). To find 



the general solution, recall that given two solutions of this second order ODE, the 
Wronskian = /1/2 — /1/2 is a constant. It follows that the general solution can 
be written in the form 

y\u) = ^i/iW/^tI^- (7-35) 
Imposing the relevant boundary conditions, we deduce 



Auiu,u') = h{u)h{u') / (7.36) 



In the same way, a particular solution for z = 2 is given from (|7.24|) as y"^ = gi{u) = 
r{u) sin6'(M). in this case, it is easy to see by inspection (in fact, as a consequence of 
angular momentum conservation) that the general solution can be written simply in 
the form 

y^{u) = cir{u) sm{e{u) + C2) . (7.37) 

This direction is obviously the focusing one since y'^{u) oscillates with a amplitude 
that decreases as the singularity is approached. The other direction de-focuses and 
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y^{u) increases as the singularity is approached. Applying the boundary conditions, 
we find 

A22iu,u') = -r{u)r{u')sm{e{u)-e{u')) . (7.38) 
This determines the VVM matrix A{u,u'). The non- vanishing components are 

du 1-1 



Aii(m,m') = iu — u 



r(u)r(u)r(u)r(u 



(7.39) 



and 



r / duM-^ 

A22{u,u') = h{u-u') r{u)r{u')sm[h / — ) • (7.40) 

L ^ Ju' ''^ 



With these expressions for the VVM matrix, the refractive index for scalar and 
spinor QED can now be evaluated by substituting into the formulae ( [4.10| ), ( ^.2^ ) 
using the explicit solution r{u) of the geodesic equation. The solution, however, 
does not show any special algebraic simplicity and an explicit numerical evaluation 
is not expected to reveal any distinctive new physics beyond what can be deduced 
by considering the particular limits and special cases discussed above. 

Rather more enlightening is the study of the refractive index and its analytic 
properties in the near-singularity limit. As shown above, the Penrose limit in this 
region is a Ricci fiat singular homogeneous plane wave with hn = — /i22 = 
The VVM matrix components are therefore given by substituting the coefficients 
«! = I and 02 = I into the general formula ( |6.13|) . This gives 



A / i\ 7{uu'y^^{u — u') A / n U — u' /^, 

^"(«'«) = 5(«^/^ - u'V^) ' ^^^("'"^ = 5(W)^/^(nV5-,./5) • (7-41) 



For scalar QED, the refractive index is given by ( |6.15| ) with the other polarization 



obtained by swapping the powers | ^ |- Notice that the affine parameter u is 
negative, with the singularity at u = 0. In this case, the integrand has a branch 
point at t = — i.e. when u' lies at the singularity, and double poles at 

t = \u\{e^'''-Pl'^ - I) , p=l,2,...,6. (7.42) 

as shown in Fig.^. 

The refractive index in the near singularity limit can be evaluated numerically 
and the results are shown in Fig.^. As expected for a Ricci fiat spacetime, the 
signs of the corrections are the same for scalar and spinor QED, with both showing 
birefringence. In this example, however, unlike the case of the Ricci fiat symmetric 
plane wave, the sign of \m.n{(jj) differs for the different polarization states. 
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Figure 4: Analytic structure of the t integrand of T^'^'^^^'^ {u; z) in the near-singularity region of 
Schwarzschild spacetime. The branch point lies at the position of the singularity and the other 
points arc double poles. 




Figure 5: n{uj) — 1 for the near singularity region of the Schwarzschild black hole plotted as a 
function of log lo: continuous (real part, polarization i ^ 1); big dashes (imaginary part, i = 1); 
small dashes (real part, polarization i = 2) and dots (imaginary part, i = 2), for scalar QED (left) 
and QED (right). 

7.3 Kerr black hole 

The case of the Kerr spacetime describing a rotating black hole can be analyzed in 
the same way and provides a very nice illustration of the simplicity of the Newman- 
Penrose tetrad approach. The Kerr metric is: 

ds^ = -p^^df + p^^dr^ + p^de^ + ^T?siYi^e{d<p-ujdtf , (7.43) 

where A = — 2Mr + a?, = (r^ + a^)^ — sin^ 9 A, p"^ = + cos^ 9 and 
u) = 2Mra/T?. The metric is specified by two parameters, M and a, where M is the 
mass and Ma the angular momentum (about the rotation axis ^ = 0) as measured 
from infinity. Note that the time component simplifies, gu = (1 " 2Mr / p^). 

The condition A(r) = 0, for which there is a coordinate singularity similar to 
that in the Schwarzschild metric, has two real solutions, r = r-t = M ± a/ M ^ — 



(assuming a < M). The larger, r = r+, is the event horizon. The region r < r_ 
contains a ring singularity. However, for the Kerr metric, this does not coincide with 
the condition = 0, which defines the boundary of the region, the ergosphere, where 
an asymptotically timelike Killing vector becomes spacelike. Within the ergosphere, 
even null curves are pulled round in the direction of the rotation. Its outer limit 
rE{0) = M + a/M^ — cos^ 6, known as the stationary limit surface, coincides with 
the event horizon at the poles and is equal to the Schwarzschild radius 2M at the 
equator. 

The null geodesies in the Kerr metric have been extensively studied (see ref. [11] 
for a comprehensive discussion, and [24] for a first analysis of photon propagation 
in QED). Although the equatorial orbits are of special interest and simplicity, our 
formalism allows us to consider immediately the general case of non-planar orbits. 
In this case, analogous to ( [7.16| ), we have A;'^ = E~^{i, r, 9, 0) where 



E / ^ . 2Mra 



a sin^ 9 — 6)^ 



^ (^r^ + (a^ -b^- qy + 2Mr ((a - hf + q) - a'g) ' 



E 

7 

E /. 2Mr 



[q + (a^ sin^ 9-1?) cot^ i 



(7.44) 



The orbits are characterized by the energy E and two impact parameters, the familiar 
h = Lz/E and q, its analogue in the ^-plane [11]. Both are conserved quantities. 

Like Schwarzschild, the Kerr spacetime is Petrov type D and the only non- 
vanishing component of the Weyl tensor is = Cemnm referred to the null basis^^ 

r = i(r2 + a^A,0,a) 

= ^(r' + a',-A,0,a) (7.45) 



M _ _! , 



m'' = (ia sin 9, 0, 1, — — -) 



\/2p ' ' ' sin 6^ 

where p = r + iacos9. In this basis, = —M/{p*Y 



16 



The equivalent covariant vectors required for the calculations here are 

£^= {-1,^,0, asinH) , n^ = ^(_i,_:^,o,asin2 
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At this point, we can simply invoke the general discussion of type D spacetimes 
in section 7.1 to write down the Penrose hmit. The Kerr spacetime is Ricci fiat, 
so from ( |7.14| ) we see that the eigenvalues of the plane wave profile function 



corresponding to the null geodesic with tangent vector k'^ are simply 

hij = 6ij{-iy\'^2\ k-ik-n, (7.46) 



with |^E'2| = M/p^. A straightforward algebraic exercise now shows that 

1 

consistent with the general theorem that 



k-ek-n = _((a-6)2 + g) , (7.47) 



\K,\^ = 41^^21'^ k-ek-n = 2M-^{a-bf + q) (7.48) 
is conserved along the geodesic. For the general Kerr orbit, we therefore have 

hu = - ^ ((« - br + , ^22 = ^ ((a - bf + q) (7.49) 

This is a remarkable simplification and illustrates very clearly the dependence of the 
Penrose limit on both the spacetime curvature and the characteristics of the chosen 
null geodesic. 

Eq.( |7.49| ) is a straightforward generalization of the Schwarzschild result ( [7.25| ) 
and similar consequences follow here. Clearly, the Penrose limit is fiat for a principal 
null geodesic, i.e. k^ = i^. Once again, there is an unstable closed orbit for a constant 
value Tc of the radial coordinate. Imposing r = 0, and specializing for simplicity to 
equatorial orbits, determines the relation Tc = 3M(6c — a) /{be + a), though here 
the expression for the critical radius is more complicated, involving the solution of a 



cubic equation, and we have Vc = 2M 



for the retrograde and 



l + cos(^|cos"^(±-5:) 

direct orbits respectively [11]. Evaluating hij for this orbit, we find the non- vanishing 
elements 

1 fbc + a\^ , 3 1 /6c + a^2 

ill ^ ^ 



^e - 'ufAb:^) ' - - ufAb:^) ^^-^^^ 

Since this is a constant, we confirm that, just as in the Schwarzschild case, the 
Penrose limit for the critical orbit in Kerr spacetime is a Ricci fiat symmetric plane 
wave. 

The near-singularity limit follows straightforwardly. Restricting to equatorial 
orbits, we determine the small r behaviour of the radial geodesic equation from 
(OD as 

= \s ' > (7-51) 
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and integrating, we find = ^-Mip — a)^u^ and therefore tfie Penrose limit profile 
function is a homogeneous plane wave identical to the Schwarzschild result (|7.33| ). 
This reflects the equivalent power-law nature of the singularities for Schwarzschild 
and Kerr black holes. This result implies that the singularity behaviour of the re- 
fractive index is the same as for the Schwarzschild solution. 



7.4 Reissner-Nordstrom black hole 



It is very simple to generalize the calculation of the Penrose Limit for the Schwarzschild 
black hole to the Reissner-Nordstrom black hole. Here we shall simply state the result 
that 

3Mb^ Ab^Q^ , 3M62 

+ —T- ' ^22 = — T -T- ■ (7-52) 



til 



^5 ^6 ' 



Notice that the hole is not Ricci fiat so that 



Ruu = hu + /i22 = • (7.53) 



This result is not as it stands useful for the calculation of the refractive index 
because in this case there is background electro-magnetic field that would have to be 
taken into account in addition to the curvature of spacetime [25] . 



8. Cosmological spacetimes 



We now consider a cosmological spacetime with an initial singularity, the Friedmann- 
Robertson- Walker universe. This example is of more than just conceptual interest 
since modifications to the speed of light in the very early universe are potentially 
important in the context of the horizon problem, one of the original motivations for 
inflation. 



The metric for a FRW spacetime has the form 



•1) 



where /^(r) = r, sinr and sinhr for k = 0, +1 and —1 respectively, representing a 
flat, closed or open universe, and dfl"^ is the usual metric for the two-sphere. With 
no loss of generality, we consider geodesies with zero angular momentum on this 
transverse space. These null geodesies are well known. We have 



. 1 

t = - 



1 
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Note that we use the notation t = ^ for -u-derivativcs; in what follows we denote 
t-derivatives by a'(t) = ^, etc. Here, u is as usual the affine parameter along the null 
geodesic 7. We therefore find = (a~^, a~^, 0, O) in a {t, r, 9, 0) coordinate basis. 

The first step in finding the Penrose hmit by the conventional method [12] de- 
scribed earlier is to identify a parallel transported frame Ea along 7. Here, we have: 

E^f^ = = (a-\a-2, 0,0) £;/ = -^(-a,l,0,0) 

^1^ = (a/,)-i(0, 0,1,0) E2'' = (a/,)-i(0,0,0,l/sin^) 

The matrix Qij = Ei^Ej^V Jzy = 1,2) follows straightforwardly and we find 
after a short calculation that^^ 

% = Sij^\og{af^) , (8.4) 

with a = a{t{uj), = /«;(r(M)) along 7. The profile function for the Penrose plane 
wave is then 

hij = -Sij{af^) ^^(«/«) • (8-5) 

This confirms hu = /i22, as expected since the Penrose limit inherits the conformally 
flat property of the original FRW metric. Then, since the metric function /^(r) 
satisfies + nf^ = for all k, we readily find [12] 



hij — S-ij 



- 4) , (8.6) 



Finally, to make contact with conventional FRW dynamics, we can re-express this in 
terms of t-derivatives, giving 



(Jj 

using the standard Friedmann and acceleration equations, 

a'2 K SttG 
-P 



(8.7) 



3 

a" 47rG 

^(P + 3p) 



(8.8) 



a 



17. 



The calculations in this section involve the following ChristofFel symbols: 

r*^ = aa' Tl^ = a-^a' 

Tle = aa'fl = aa' flsm^ 6 ^6 = -/% ^U = -f%^^^ 
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and the equation of state p = wp. 



The Newman-Penrose null tetrad derivation of the Penrose limit is especially 
quick in this example. Recall from the discussion in the preceding section that the 
profile function is given by (|7.8| ), where here we can identify the NP basis vectors 
L, N, M, M with the of (|8.3|) . Since the FRW spacetime is conformally flat, we 
immediately have 



To check this with the previous result ( |8.7| ), recall that the Ricci tensor is given 
by Rfj,u = SnGiT^u — \T^\g^v) where the FRW energy-momentum tensor is T^^ = 
{p + p)t^ty + pg^y with a unit vector in the t-direction. It follows immediately that 

= %TTG{p + p){Eu-tf (8.10) 
= dij^^{p + p) , 

since the terms proportional to g^v do not contribute as Eu'^ is null. 

To take this further, we now use the standard conservation equation 

^ipa') = -3wpa^ , (8.11) 
da 

which follows from V ^T'^'^ = 0, to deduce as usual p ~ a~'^^^'^\ It is convenient 
here to introduce the parameter 7 = ^-^^^j ■ The three standard equations of state 
considered in FRW cosmology are = (7 = 2/3) for non-relativistic matter, 
w = 1/3 (7 = 1/2) for radiation, and w = — 1 (7 — >■ 00) for a cosmological constant. 
The latter implies a fiat Penrose limit, as is evident from (|8.7|) , since a vacuum 
spacetime with a cosmological constant, i.e. de Sitter space, is maximally symmet- 
ric. Combining ( 8.11|) with the Friedmann equations, and introducing the constant 



C-y^y^poT, we now have 

, 2 (8.12) 

a" = C,(l-i)aM 

showing clearly that 7 < 1 (w < — 1/3) is the threshold for an accelerating universe. 
We therefore find from (|8.7| ) that the Penrose limit profile function is 

hij = SijC^-f-^a'^^^^^^ , (8.13) 
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where the dependence of a{u) on u follows from ( ^.12 ) as 



u = / da , . (8.14) 

2(1-1) 

Cya ^ -y' — K 



It is now simplest to consider the various possible cases separately, 
(i) Spatially flat universe, k = 0, for any 7: 



In this case, the integral in (|8.14| ) is trivial and the dependence of a{u) on u is 

7+1 



given implicitly by u = C-y ^^'^ -^-^■a^'^'^ \ The profile function ( p.l3|) follows immedi 



ately: 



Perhaps unsurprisingly given our experience with black hole spacetimes with singu- 
larities in the previous section, we find here that the Penrose limit of a flat FRW 
spacetime is a conformally fiat singular homogeneous plane wave with coefficient 

(ii) 7 = 1, for any k: 

As already noted, an equation of state with 7 = 1 (ly = — |) is on the boundary 
between a decelerating (0 < 7 < 1) and accelerating (7 > 1) universe. From the 
Friedmann equation, we see this requires Ci — k > 0, which only becomes a constraint 
on the energy density for a closed universe, k = 1. In this case, we find u = | 
and the profile function can be expressed as 

Again we find a singular homogeneous plane wave, though with a K-dependent mod- 



ification of the coefficient in ( |8.15| ). The case k = 1 is interesting in that in this case 
a^ = (l — Ci)~^<0 and so at is imaginary and, following the discussion in Section 
O, there are an infinite sequence of conjugate points. 



In these special cases, the characteristic 1/u^ behaviour of hij holds exactly for 
all u, i.e. along the whole geodesic. The following results, in contrast, describe only 
the leading behaviour in the near-singularity limit. 

(iii) Decelerating universe, < 7 < 1, for all k: 
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Here, the first term in the denominator of (^.14 ) dominates for small a, and the 



dependence of a{u) on u is exactly as in case(i). So again we find, for all k, 

7 1 

in the near-singularity limit. Notice that for the case of a matter dominated universe, 
7 = I (ti; = 0), the coefficient in ( ^.17] ) is ^, precisely the same as for the near- 



singularity limit of Schwarzschild, or Kerr, spacetime [12]. 

(iv) Accelerating universe, 7 > 1: 

For an accelerating, 7 > 1, universe, we have already found the exact solution 



( p.l5| ) for K = 0. For k = 1, the Friedmann equation does not have a solution if a 
becomes too small - the parameters k = 1,7 > 1 do not correspond to a spacetime 
with an initial singularity. This leaves the open universe k = — 1. This time, for 



sufficiently small a, the second term in the denominator of (|8.14|) dominates and we 



find u ~ ^a^, leading to 



h,, = <5,^.2-(i+^)7-iC^-^ . (8.1^ 

u -> 



The singular behaviour of hij is softer than 1/u in this case; this can be interpreted 
as due to a cancellation of the leading singular behaviour from the two terms 
in ( ^.61) when u = \a^. 



In order to consider the refractive index in these cosmological examples, the 
following conceptual issue arises. Notice that the expression for the index in ( [4.1^ ) 
and (|4.14|) , depends on the past history of the geodesic. The intuitive understanding 
is that a definition of the refractive index involves waves coming in from past null 
infinity. In an FRW universe, the singularity cuts off this past history and the 
refractive index cannot, as such, be defined in this way. However, another way to 
think about the vacuum polarization is in terms of an initial value problem: rather 
than have the waves come in from past null infinity, one defines a wave at some initial 
time. In this case one should think in terms of the retarded Green function rather 
than the refractive index. This quantity was considered in [1]. Effectively defining 
an initial value problem provides an the upper limit for the t integral in ( [4.14| ) of the 
form mq, where the initial value surface is at m = Uq. 

'•"0-^^ dt 







1 - A(u,u -t)J det A (m, u - t) 



119) 



The resulting quantity n{u; 00) then plays the role of an effective refractive index for 
waves that begin on the initial value surface which, therefore, depends implicitly on 
the choice of the initial value surface. 
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The presence of the cut-off on the t integral changes the low frequency behaviour 
of the refractive index from the expressions in (|2.17|) and ( p. 181) to 



n-f -(n, no; = ^^. + ^^^(^^n(^)<^,, + 2/?..„,(t.))F( ^^^^^_ ^^^ ) +■ ■ ■ , (8.20) 

while for spinor QED, 

n^'''"\u,Uo]uj) = 6ij - {l?>Ruu{,u)5ij - 4R^iuj{u))F( ^ ) H , 

■' 1807rm^ \m^{u — Un)'' 



where ^ 

F{x) = 1 - 6 / - Oe-*/('^('-^)") . 

Jo 

The effect is that the low-frequency behaviour oscillates as a function of uj. 



(8.21) 
(8.22) 



As an example of the kind of behaviour one encounters with this definition of 
the refractive index, consider the Milne universe [26] which has spatially flat sections 
K = and metric 

ds^ = -dt^ + aH^idx^ + dy^ + dz^) . (8.23) 

A congruence of null geodesies is simply obtained by taking y and z constant along 
with t = atx, which is solved by taking x = (2a)~^log'U — Xq, for constant Xq. The 
variables {t, xq, y, z) form a set of Rosen coordinates for the congruence. Writing the 
metric in terms of these coordinates we have 

ds^ = -2at dt dxo + aH'^{dxl + dy"^ + dz^) . (8.24) 

In order to take the Penrose limit, we choose an affine parameter u = at'^/2 and 
identify xq = V , = y and = z. Then the Penrose is easily taken giving, in 
both Rosen and Brinkmann coordinates, 

ds^ = -2du dV + 2au {dY^ dY^ + dY^ dY^) 

1 (8.25) 
= -2dudv - -^((y'? + {y^f)du^ + dy'dy' . 

This is precisely a singular homogeneous plane wave of the type discussed in Section 
6^ with aj = and the VVM matrix may be extracted from (|6.13|) by taking the 



limit — i> 0: 

A,,(n,«') = 5,, ^'"^ . (8.26) 

Numerical evaluations of the refractive index are shown in figure ^ for scalar and 
spinor QED. In this case the spacetime is conformally flat and so there is no bi- 
refringence. At low frequency we see the characteristic oscillations that match ( p.20| ) 
and {^^. 
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Figure 6: Ren — 1 (continuous) and Imn (dashed) plotted as a function of logw in suitable units 
for the Milne universe for (a) scalar QED and (b) spinor QED. 



9. Gravitational Waves 



Finally, we consider another major class of spacetimes - gravitational waves. Two 
examples are studied in detail, viz. weak gravitational waves and gravitational shock- 
waves, the latter being of particular interest because of the use of the Aichelburg-Sexl 
metric [27] in investigations of Planck energy scattering in quantum gravity. 



9.1 Weak gravitational waves 



The spacetime metric for a weak gravitational wave has the following form in Rosen 
coordinates: 

ds'^ = - 2du dV + {1 + e cos uu)dY^ dY^ + [l-e cos uujdY^dY^ . (9.1) 

Here, and in the following, e is small and we work to linear order. The transformation 
to Brinkmann coordinates is made using the zweibein 

Ea'iu)=6:{l-i-iy^cosuu) , (9.2) 

giving 

ds^ = -2dudv - —cos uu[{yY ~ {yY]du^ + dy^ dy^ . (9.3) 

Since hn = — /i22 = |ez/^cosz/M, this is a Ricci flat spacetime and, as usual, will 
display gravitational birefringence. 



The metric (|9.3| ) is already in the form of the Penrose limit for a photon travelling 
in the opposite direction to the gravitational wave along the null geodesic given by 
( |3.30| ). Notice that a photon travelling in the same direction as the gravitational wave 



has a flat Penrose limit and there is no vacuum polarization effect on its refractive 
index [4,28]. 
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The Jacobi fields therefore satisfy 

^2 (±) ,,.2 

^±^cosz.«yW = 0, (9.4) 

where we have labelled the two polarizations y^^^ = and y^~^ = y'^. This can 
easily be solved perturbatively in e, with solution to linear order 

y^'^^u) = Ci + UC2 ± -(ci + UC2) cos uu =F — - sin i/u . (9.5) 

The (diagonal) matrix A{u,u') is the solution of the Jacobi equation subject to the 
now-familiar boundary conditions and we find its components for the tow polariza- 
tions to be 

^{±){u, u ) = u — u 



± ^ ( cos vu + cos z/u') T — ( sin uu — sin z/m') . 



This determines the eigenvalues of the Van-Vleck Morette matrix: 

-(cos z/u + cos z/u') ± — 

2^ ' v{u-u') 



A(j-)(u,u') = 1 =1= -(cos z/u + cos z/u') ± — j^(sin z/M — sin z/m') . (9.7) 



Notice that this is regular for all m, m', reflecting the fact that in this weak field (small 
e) approximation there are no conjugate points in the null geodesic congruence. 



The refractive index for scalar QED is given as usual by ([4.10|) where for the 2 
polarizations we have 







X 

=Fe \_f{z) cos uu + g{z) sin z/-u] 



-(cos z/M -|- cos z/('u — t)) (sin uu — sin uiu — t)) 

2 ut 



{9.t 



with 



f{z) = ^ [2z/(l + log z) + {z- v) \og{z - z/) - (2; + z/) log(2; + v)\ , 
(7(2;) = ^ ['^^ + 2-2^ log z -|- z{y — z) \og{z — u) — z{z + u) \og{z + z/)] 



(9.9) 



For spinor QED the only difference relative to scalar QED is an overall multiple of 
4.18 

J^l^^^^\u;z)=ATl^;-^iu;z) , (9.10) 



-'^^In fact this generalizes to any Ricci flat background: the spinor and scalar QED results to linear 
order in the curvature, but allowing any number of derivatives, only differ by a factor of 4. 
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and so we will continue to discuss the scalar case. 



The functions f{z), g{z) have branch points at 0, oo and z = ±z/ and this 
means that n(^±){u;uj) will have branch points at ±00 and ±2m^/z/. In particular, 
the branch points at ±2m^/z/ are points of non-analyticity of the refractive index. 
This non-analyticity manifests itself by the fact that lmf{z) and lie g{z) are zero 
for 2; G R > z/, while for 2; G R < z/, 

lmf{z) = Reg{z) = ^^^^^^ . (9.11) 



The low frequency expansion of the refractive index follows immediately from 
and we find: 



n^'^\u; uj) 



1 ± 



1 1 /LUU 

— + I — 



L360 6300 Vm2 
.aez/2 r 1 (^^\ 1 i'uju\^ 
40 W/ ^ 10395 W 



while at high frequencies, 



1 ± t sm uu + 



cos uu 



sm uu 



(9.12) 



(9.13) 



Notice that the low frequency expansion is not sensitive to the non-perturbative 
contribution ( 9.11| ). The leading terms for both real and imaginary parts can however 
be read off from the effective Lagrangian described in section 2. 





Figure 7: (a) Ren — 1 and (b) Iiiiti for u = 0.2, m = ly = 1 plotted as a function of log vui/m? in 
units of ea for both polarizations. The point of non-analyticity at w = ^rr? jv is quite clear. (Note 
that the fact that the polarizations do not quite give mirror images is an artifact of the numerical 
approximation. 



The full form of the frequency dependence of the real and imaginary parts of 
the refractive index is plotted numerically in Fig. |^, evaluated at a fixed point on 
the photon trajectory, for both scalar and spinor QED. For the first polarization. 
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has a Ren(co') which quahtatively looks similar to a conventional dielectric with a 
single characteristic frequency uj ~ m^/z/, together with an imaginary part Imn(ci;) 
which is negative in the perturbative small uj regime before being dominated by the 
positive non-perturbative contribution ( 9.11| ) above the branch point at = 2w? jv. 
Since this is a Ricci flat spacetime, the second polarization has the opposite sign, 
being superluminal at low frequencies and with \m.n[LS) < above the branch point, 
normally an indication of amplification rather than dispersive scattering. The roles 
of the two polarizations of course change along the photon trajectory through the 
background gravitational wave. 



9.2 Gravitational Shockwaves 



A gravitational shockwave is described by the Aichelburg-Sexl metric [27] 

ds^ = ~2du dv + f{r)5{u) du^ + dx^ + dy"^ , (9.14) 

where x = rcos0 and y = rsin0. This describes an axis symmetric, plane-fronted 
wave moving at the speed of light, with a profile function /(r) in the transverse 
direction determined by the nature of the source: a matter source Tuu = pij)5{u) 
corresponds to a profile function satisfying A/(r) = — 167rp(r), where A is the 2-dim 
Laplacian. We consider two sources of special interest: a particle with p(r) = fiS{x), 
which gives a profile function /(r) = — 4/ilogr^, and a homogeneous beam p(r) = /i, 
for which /(r) = —Anfir^. 

The AS metric has many applications in general relativity and more recently in 
quantum gravity and string theory, where it has been widely used in investigations 
of Planck energy scattering (see refs. [29-31] for a selection of papers). Recently, 
string scattering in a classical AS metric has been studied [32] as a model of Planck 
scattering, and it would be interesting to compare these results, characterized by 
the string scale, with the equivalent QFT scattering, where the scale is set by the 
Compton wavelength of the particle in the vacuum polarization loop. For now, we 
focus on the effect of a shockwave collision on the refractive index in QED. 

The null geodesies describing a photon travelling in the opposite direction to the 
gravitational shockwave are well-known; those corresponding to a homogeneous beam 
were first discussed by Ferrari, Pendenza and Veneziano in ref. [33] (see also [34]). 
With u as the affine parameter, they satisfy 

r = R+\f\R)ue{u) , 

v = V + \f{R)e{u) + \f\Rfue{u) , (9.15) 

= $, 
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where (V, R, $) are constants which label a congruence of null geodesies. Conse- 
quently, these variables are a natural set of Rosen coordinates. It is now straightfor- 
ward to change variables from (m, v, r, (p) to (m, V, R, $): 



dr = dR+ \f"{R)u9(u)dR + \f\R)9{u)du, 



dv = dV + \f\R)e{u)dR + \f{R)5{u)du 

+ \f\R)f"{R)ue{u)dR+\f{Rfe{u)du 



d(f) = d^ , 

where we used the consistent replacement u6{u) — * 0. The metric in Rosen coordi- 
nates is therefore 

ds^ = -2dudV+ (l + lf"{R)ue{u)ydR^ + [R+\f\R)ue{u)yd^^ . (9.17) 



We can now calculate the Penrose limit around a geodesic with R = Rq hy 
defining shifted Rosen coordinates 

Y^ = R-Ro, F^ = i?o$, (9.18) 

The metric is now in standard form and it is trivial to implement the Penrose limit; 
we find 

df = -2dudV + (l + y\Ro)u9{u)y{dY^f + + lRo^f'{Ro)ue{u)y{dY^f . 

(9.19) 

The zweibein E^a{u) can be read off directly (see ( p.27| )), and has non- vanishing 
components 

E\{u) = 1 + \f"{R^)u9{u) , E\{u) = 1 + \R^^f{R^)u9{u) . (9.20) 

We therefore find the Penrose limit in Brinkmann coordinates, with a diagonal profile 
function hij given by ( p.31| ) as 

hiiu) = -y"{R,)S{u) , h22{u) = -lR^'f'{Ro)S{u) . (9.21) 

For the "particle" AS metric with /(r) = — 4/xlogr^ we have 

h.,{u) = {-iy5.,%5{u) , (9.22) 

which is an example of a Ricci flat plane wave. The related "beam" AS metric with 
/(r) = — yur^ gives 

hij{u) = 6ijfi6{u) , (9.23) 
which is a conformally flat shockwave. 
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The Van Vleck-Morette matrix is found as usual by solving the Jacobi equations 
subject to the appropriate boundary conditions, with hij given by (|9.22| ), (|9.23| ). 
Writing the above zweibeins as E^a = (l ~ \u6{u))5'^a^ so that h^j = Xi^ij, we find, 
in general, if u and u' are on different sides of the shockwave, uu' < 0, 



\u — u" 



^iji^,^') = I — hr—i — Aj^ (9-24) 

\u — u'\ + Xiuu' 

while Aij{u,u') = 1 if u,u' are on the same side of the shockwave, uu' > 0. When 
uu' < for a polarization with Aj > 0, there are conjugate points with |m| + \u'\ = 
Xi\uu'\ provided the gravitational field of the shockwave is sufficiently strong, Xi\u\ > 
1 and Ajl^'l > 1. In fact, the relation for conjugate points is 

— + — = Xi, 9.25 

m p I 

which is the lens formula for a lens of focal length A^"^. When the Aj > the lens 
is converging and a conjugate point exists as long as u and u' are large enough, as 
above. In this case, the refractive index integral will have a branch point on the real 
t axis. Correspondingly, when Aj < the lens is diverging and no conjugate points 
exist. 

The refractive index is given by the formulae in sections 4 ad 5 and is easily 
evaluated for both the particle (Ricci-fiat) and beam (conformally fiat) AS space- 
times with either scalar or spinor QED. For example, consider the conformally fiat 
shockwave and scalar QED. We have, for either polarization. 



Hu-z) = e{u) / -^te-'^' 1- ( 



t 



-t + Xu{u — t) 



e{u) [ - zr(-l, ^uz) + - 1, ^) 



(9.26) 



where T{a,z) = f^'^e'^ dt is the incomplete gamma function. Notice that the 
lower limit of the integral is at t = u, in order that the point m' = -u — t < is on 
the other side of the shockwave. The factor of 6{u) ensures that the shockwave lies 
in the past of the point u. 

In this case, A is negative and the integrand of J-'{u; z) has no poles or branch 
points on the real t axis. Nevertheless, inspection of the explicit solution in terms of 
gamma functions makes it clear that J^{u\ z) itself still has branch cuts from z = to 
oo. The refractive index n{u) evaluated numerically from ( p.26|) is shown in Fig. ^ 



The small u behaviour of the refractive index can be obtained analytically by 
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Figure 8: n — 1 (real part continuous, imaginary part dashed) for A = 1 and u — 2 (left) and 
u — ^ (right) plotted as a function of logo; for scalar QED. 



expanding the integrand of J-'+{u; z) in powers oit — u. In all cases, 

(9.27) 



J u 

= i9{u)aoe-''^{izuy^-^T{l + 7) + ■ ■ ■ , 



where, depending on the case, 7 = or 1 and is some Am dependent constant. 
The ^ integral is then dominated by a stationary phase with ^ = ^ which gives a 
behaviour 

niu;uj) = l-aC(M,m,A)^(M)cuT+^/V2^™'"/" , (9.28) 

where C{u,m, A) is some uj independent factor. This shows clearly how in this case 
the non-vanishing imaginary part of the refractive index arises from the oscillating 
phase factor in the integrand evaluated at the lower limit of the t-integral, even 
though there are no singularities or branch points off the real axis. We can make 
these formulae even more explicit for the conformally flat scalar QED case. In this 
case from ( |9.26|) , we get 



On the other hand, for large u, 

n{u-u) = I- 9{u)+\ogicco)o(^) . (9.30) 

Note the oscillations at small u as shown in the numerical evaluation in Fig. |. 

It is worth pointing out that in this case we cannot match the low frequency 
behaviour to the effective action because the curvature is not differentiable. The other 
thing this computation makes clear is the non-local natural of vacuum polarization: 
the refractive index is non-trivial even in a region of fiat space because the e"*"e" loop 
can extend far enough back in time to touch the shockwave. 
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10. Discussion 



The propagation of light, and other quantum fields, is a fundamental and difficult 
problem for QFT in curved spacetime. In particular, it has raised the question of how 
the original Drummond and Hathrell discovery of low-frequency superluminal prop- 
agation could be consistent with causality. In our previous work, we have resolved 
this issue by developing techniques which allow an explicit computation of the full 
frequency dependence of the refractive index and Green functions for (scalar) QED 
in curved spacetime. 

In the present paper, wc have further developed the "phenomenology" of the re- 
fractive index in a number of important spacetimes, namely black holes, gravitational 
waves and FRW universes, and extended our analysis to spinor QED. A crucial ob- 
servation is that, to leading order in R/rri^, the refractive index depends only on the 
Penrose limit of the background spacetime around the original null geodesic. This is 
a very powerful simplification and reveals a number of universal features. For exam- 
ple, propagation in many examples reduces in the Penrose limit to the simple case of 
homogeneous plane waves, especially the singular homogeneous plane waves which 
are associated with the near-singularity limits for both black holes and cosmological 
spacetimes. 

A key element in understanding the causal properties of light propagation, and 
more generally the retarded and advanced Green functions, is the analytic structure 
of the refractive index. As we have shown, this is in turn controlled by the Van Vleck- 
Morette matrix which describes the geometry of geodesic deviation. The examples 
we have studied here display a rich variety of analytic structures, very different 
from the simple picture in flat spacetime. The Green functions and refractive index 
therefore display many interesting and physically important features special to curved 
spacetimes; in particular, conventional dispersion relation theory in QFT, exemplified 
by the Kramers- Kronig relation, is radically changed. Nevertheless, in all cases we 
have investigated, the analytic structure still allows the definition of retarded and 
advanced Green functions with the required causal properties. We can conclude that 
QED in curved spacetime is indeed a perfectly causal theory despite the unusual and 
remarkable occurrence of superluminal phase velocities at low frequencies. 

Whilst this goal has been achieved and one can finally see explicitly how a low- 
frequency phase velocity bigger than c is compatible with retarded Green functions 
that vanish outside the lightcone, a new puzzle arises in that the imaginary part 
of the refractive index is found in many examples to be negative. This indicates 
that the optical theorem is violated in curved spacetime. A full explanation of this 
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phenomenon will be presented elsewhere [7]. 
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